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ABSTRACT 

Spatial extension is an important characteristic for correctly associating 7-ray- 
emitting sources with their counterparts at other wavelengths and for obtaining 
an unbiased model of their spectra. We present a new method for quantifying 
the spatial extension of sources detected by the Large Area Telescope (LAT), the 
primary science instrument on the Fermi Gamma-ray Space Telescope {Fermi). 
We perform a series of Monte Carlo simulations to validate this tool and calculate 
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the LAT threshold for detecting the spatial extension of sources. We then test 
all sources in the second Fermi-LAT catalog (2FGL) for extension. We report 
the detection of seven new spatially extended sources. 

s Subject headings: Catalogs; Fermi Gamma-ray Space Telescope; Gamma rays: 

9 observations; ISM: supernova remnants; Methods: statistical; pulsar wind nebula 
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A number of astrophysical source classes including supernova remnants (SNRs), pulsar 
wind nebulae (PWNe), molecular clouds, normal galaxies, and galaxy clusters are expected 
to be spatially resolvable by the Large Area Telescope (LAT), the primary instrument on 
the Fermi Gamma-ray Space Telescope (Fe rmi). Additionally, da rk matter satellites are also 
hypothesized to be spatially extended. See lAtwood et al.l (|2009a ) for pre-launch predictions. 
The LAT has detected seven SNRs which are significantly extende d at GeV energies: W51C 



W30, I C 4 43, W28, W44, RX J1713.7 -3946, and the Cygnus Loop flAbdo et al.l2009all2012b 



2010hffl l201ll ; Katagiri et al.ll201lf ). In addition, three extended PWN ha ve been detected 



by th e LAT: MSH 15-52, Vela X, and HESS J1825- 137 f lAbdo et al.ll2010aM iGrondin et al. 



201ll ). Two nearby galaxies, the Large and Small Magellanic Cloud s, and the lobes o f one 
radio galaxy, Centaurus A, were spatially resolved at GeV energies ( Abdo et al. 2010i b.cl) . 
A number of additional sources detected at GeV energies are positionally coincident with 
sources that exhibit large enough extension at other wavelengths to be spatially resolvable 
by the LAT at GeV energies. In particular, there are 59 GeV sources in the s econd Fermi 



Sourc e Catalog (2FGL) that might be associated with extended SNRs (2FGL, iNolan et al. 



20121 ). Previous analyses of extended LAT sources were performed as dedicated studies of 
individual sources so we expect that a systematic scan of all LAT-detected sources could 
uncover additional spatially extended sources. 

The current generation of air Cherenkov detectors have made it apparent that many 
sources can be spatially resolved at even higher energies. Most prominent was a survey of 
the Galactic plane using the High Energy Stereoscopic System (H.E.S.S) which reported 14 
spati ally extended sources with extensions varying from ~ 0?1 to ~ 0?25 ( lAharonian et al. 



20061 ). Within ou r Galaxy very few source s detected at TeV energies, most notably the 7-ray 



binar ies LS 5039 (lAharonian et al.ll2006al) LS 1 +61-303 (Albert et al. 



20061: lAcciari et al. 



201 lh . HESS J0632+057 (lAharonian et al.ll2007ch . and the Crab nebula (IWeekes et al.lll989h . 
have no detectable extension. High-energy 7-rays from TeV sources are produced by the de- 
cay of 7r°s produced by hadronic interactions with interstellar matter and by relativistic 
electrons due to Inverse Compton (IC) scattering and bremsstrahlung radiation. It is plausi- 
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39 ble that the GeV and TeV emission from these sources originates from the same population 

40 of high-energy particles and so at least some of these sources should be detectable at GeV 

41 energies. Studying these TeV sources at GeV energies would help to determine the emission 

42 mechanisms producing these high energy photons. 
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The LAT is a pair conversion telescope that has been surveying the 7-ray sky since 2008 
August. The LAT has broad energy coverage (20 MeV to > 300 GeV), wide field of view 
(~ 2.4 sr), and large effective area (~ 8000 c m 2 at > 1 GeV) Addi tional information about 
the performance of the LAT can be found in lAtwood et al.l ( 12009bl ). 



Using 2 years of all-sky survey data, the LAT Collaboration published 2FGL (2FGL, 
Nolan et al.l 120121 ) . The possible counterparts of many of these sources can be spatially 
resolved when observed at other frequencies. But detecting the spatial extension of these 
sources at GeV energies is difficult because the size of the point-spread function (PSF) of 
the LAT is comparable to the typical size of many of these sources. 

The capability to spatially resolve GeV 7-ray sources is important for several reasons. 
Finding a coherent source extension across different energy bands can help to associate a LAT 
source to an otherwise confused counterpart. Furthermore, 7-ray emission from dark matter 
annihilation has been predicted to be detectable by the LAT. Some of the dark matte r 



substructure in our Galaxy could be spatially resolvable by the LAT ( iBaltz et al.ll2008l ). 
Characterization of spatial extension could help to identify this substructure. Also, due to 
the strong energy dependence of the LAT PSF, the spatial and spectral characterization of a 
source cannot be decoupled. An inaccurate spatial model will bias the spectral model of the 
source and vice versa. Specifically, modeling a spatially extended source as point-like will 
systematically soften measured spectra. Furthermore, correctly modeling source extension 
is important for understanding an entire region of the sky. For example, an imperfect 



(Abdo et al. 


2010d; 


Nolan et al. 


2012) 



measured spectra of neighboring sources in the densely populated Galactic plane. 

For these reasons, in Section [2] we present a new systematic method for analyzing spa- 
tially extended LAT sources. In Section [31 we demonstrate that this method can be used to 
test the statistical significance of the extension of a LAT source and we assess the expected 
level of bias introduced by assuming an incorrect spatial model. In Section HJ we calculate 
the LAT detection threshold to resolve the extension of a source. In Section [51 we study the 
ability of the LAT to distinguish between a single extended source and unresolved closely- 
spaced point-like sources In Section ® we further demonstrate that our detection method 
does not misidentify point-like sources as being extended by testing the extension of active 
Galactic nuclei (AGN) believed to be unresolvable. In Section [71 we systematically reanalyze 
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75 the twelve extended sources included in the 2FGL catalog and in Section [8] we describe a 

76 way to estimate systematic errors on the measured extension of a source. In Section [91 we 

77 describe a search for new spatially extended LAT sources. Finally, in Section [TU] we present 

78 the detection of the extension of nine spatially extended sources that were reported in the 

79 2FGL catalog but treated as point-like in the analysis. Two of these sources have been 
so previously analyzed in dedicated publications. 



2. Analysis Methods 



82 Morphological studies of sources using the LAT are challenging because of the strongly 

83 energy- dependent PSF that is comparable in size to the extension of many sources expected 

84 to be detected at GeV energies. Additional complications arise for sources along the Galactic 
as plane due to systematic uncertainties in the model for Galactic diffuse emission. 

se For energies below ~300 MeV, the angular resolution is limited by multiple scattering in 

87 the silicon strip tracking section of the detector and is several degrees at 100 MeV. The PSF 

as improves with energy approaching a 68% containment radius of ~ 0?2 at the highest energies 

89 (when averaged over the ac ceptance of the LAT) and is limited by the ratio of the strip pitch 

90 to the height of the tracker jAtwood et all2009bt lAbdo et al]l2009dl . l2012ah However, since 

91 most high energy astrophysical sources have spectra that decrease rapidly with increasing 

92 energy, there are typically fewer higher energy photons with improved angular resolution. 

93 Therefore sophisticated analysis techniques are required to maximize the sensitivity of the 

94 LAT to extended sources. 



2.1. Modeling Extended Sources in the pointlike Package 



A new maximum-likelihood analysis tool has been developed to address the unique 
requirements for studying spatially extended sources with the LAT. It works by maximizing 
the Poisson likelihood to detect the observed distributions of 7-rays (referred to as counts) 
given a parametrized spatial and spectral model of the sky. The da ta are binned spatially, 
using a HEALPix pixellization, and spectrally ( iGorski et al.l 120051 ) and the likelihood is 
maximized over all bins in a region. The extension of a source can be modeled by a geometric 
shape (e.g. a disk or a two-dimensional Gaussian) and the position, extension, and spectrum 



^^More information about the performance of the LAT can be found at the Fermi Science Support Center 
(FSSC, |http : //f ermi . gsf c .nasa. gov) . 
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of the source can be simultaneously fit. 

This type of analysis is unwieldy using the standard LAT likelihood analysis tool gtlike§| 
because it can only fit the spectral parameters of the model unless a more sophisticated 
iterative procedure is used. We note that gtl ike has been used in t he past in several studies 



of source extension in the LAT Collaboration ( lAbdo et al.ll2010il JbUd. l2009al ). In these studies 



a set of gtlike maximum likelihood fits at fixed extensions was used to build a profile of the 
likelihood as a function of extension. The gtlike likelihood profile approach has been shown 
to correctly reprod uce the extensio n of simulated extended sources assuming that the true 
position is known ( iGiordanol 1201 ll ). But it is not optimal because the position, extension, 
and spectrum of the source must be simultaneously fit to find the best fit parameters and 
to maximize the statistical significance of the detection. Furthermore, because the gtlike 
approach is computationally intensive, no large-scale Monte Carlo simulations have been run 
to calculate its false detection rate. 

The approach presented here is based on a second maximum likelihood fitting pa ckage 
developed in the LAT Collaboration called point like (lAbdo et al.ll2010dl ; lKerrll201ll ). The 
choice to base the spatial extension fitting on point like rather than gtlike was made due 
to considerations of computing time. The pointlike algorithm was optimized for speed 
to handle larger numbers of sources efficiently, which is important for our catalog scan 
and for being able to perform large-scale Monte Ca rlo simulatio ns to validate the analysis. 
Details on the pointlike package can be found in iKerrl (120 111 ). We extended the code to 
allow a simultaneous fit of the source extension together with the position and the spectral 
parameters. 



125 2.2. Extension Fitting 

126 In point like, one can fit the position and extension of a source under the assumption 

127 that the source model can be factorized: M(x,y,E) = S(x,y) x X(E), where S(x,y) is 

128 the spatial distribution and X(E) is the spectral distribution. To fit an extended source, 

129 pointlike convolves the ext ended source shape with the PSF (as a function of energy) 



130 and uses the minuit library (I James fc Rooslll975l ) to maximize the likelihood by simulta 



131 neously varying the position, extension, and spectrum of the source. As will be described 

132 in Section 13. 1\ simultaneously fitting the position, extension, and spectrum is important to 

133 maximize the statistical significance of the detection of the extension of a source. To avoid 

134 projection effects, the longitude and latitude of the source are not directly fit but instead 



2 gtlike is distributed publicly by the FSSC. 
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135 the displacement of the source in a reference frame centered on the source. 



The significance of the extension of a source can be calculated from the likelihood-ratio 
test. The likelihood ratio defines the test statistic (TS) by comparing the likelihood of a 
simpler hypothesis to a more complicated one: 

TS = 21og(£(if 1 )/£(#o)), (1) 

where Hi is the more complicated hypothesis and H the simpler one. For the case of the 
extension test, we compare the likelihood when assuming the source has either a point-like 
or spatially extended spatial model: 



TS cxt = 21og(£ ext /£ 



psj 



(2) 



136 

137 
138 



point like calculates TS cxt by fitting a source first with a spatially extended model and 
then as a point-like source. The interpretation of TS cx t in terms of a statistical significance 
is discussed in Section 13.11 

For extended sources with an assumed radially-symmetric shape, we optimized the cal- 
culation by performing one of the integrals analytically. The expected photon distribution 
can be written as 



PDF (r) 



PSF(|r-r'|)/ 8rc (f)r'dr / 



(3) 



where r represents the position in the sky and I STC (r) is the spatial distribution of the source. 
The PSF of the LAT is cur rently parameterize d in the Pass 7_V6 ( P7_V6) Sour ce Instrument 
Response Function (IRFs. lAbdo et al.ll2012al ) by a King function (|Kingill962l ): 



PSF(r) 



2na 2 



i-I 1+ « 

7/ V 7 



(4) 



where u = (r/a) 2 /2 and a and 7 are free parameters (jKerrll201ll ). For radially-symmetric 
extended sources, the angular part of the integral can be evaluated analytically 

»2tt 

(5) 



PDF(u) = / r'dr'I src (v) / d(f)'PSF (J 2a 2 (u + v - 2v^cos(0 - 0'))) 
Jo Jo 



dvl src (v) 



7-1 

7 



7 



7 + u + v 



7 / 1 

x 2 F 1 7/2, - 



7 



,1, 



Auv 



(7 + u + f) 5 



(6) 
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where v = (r'/a) 2 /2 and 2-^1 is the Gaussian hypergeometric function. This convolution 
formula reduces the expected photon distribution to a single numerical integral. 

There will always be a small numerical discrepancy between the expected photon dis- 
tribution derived from a true point-like source and a very small extended source due to 
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numerical error in the convolution. In most situations, this error is insignificant. But in 
particular for very bright sources, this numerical error has the potential to bias the TS for 
the extension test. Therefore, when calculating TS ext , we compare the likelihood fitting the 
source with an extended spatial model to the likelihood when the extension is fixed to a very 
small value (10~ 10 degrees in radius for a uniform disk model). 

We estimate the error on the extension of a source by fixing the position of the source 
and varying t he extension until the log of the likelihood has decreased by 1/2, corresponding 
to a la error (jEadie et al.lll97ll ). Figure [T] demonstrates this method by showing the change 
in the log of the likelihood when varying the modeled extension of the SNR IC 443. The 
localization error is calculated by fixing the extension and spectrum of the source to their 
best fit values and then fitting the log of the likelihood to a 2D Gaussian as a fun ction of 
position. This localization error algorithm is further described in iNolan et al.l (120121 ). 



155 2.3. gtlike Analysis Validation 

156 point like is important for analyses of LAT data that require many iterations such as 

157 source localization and extension fitting. On the other hand, because gtlike makes fewer 
i5B approximations in calculating the likelihood we expect the spectral parameters found with 
159 gtlike to be slightly more accurate. Furthermore, because gtlike is the standard likelihood 
wo analysis package for LAT data, it has been more extensively validated for spectral analysis, 
lei For those reasons, in the following analysis we used point like to determine the position and 

162 extension of a source and subsequently derived the spectrum using gtlike. Both gtlike and 

163 point like can be used to estimate the statistical significance of the extension of a source 

164 and we required that both methods agree for a source to be considered extended. There was 

165 good agreement between the two methods. Unless explicitly mentioned, all TS, TS ex t, and 
we spectral parameters were calculated using gtlike with the best-fit positions and extension 

167 found by pointlike. 

168 2.4. Comparing Source Sizes 

We considered two models for the surface brightness profile for extended sources: a 2D 
Gaussian model 

^GaussianO, y) = ^ exp (-(x 2 + ?/ 2 )/2cr 2 ) (7) 
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or a uniform disk model 

hi S k{x, y) 



x 2 + y 2 < a 2 



x 2 + y 2 > a 2 . 

Although these shapes are significantly different, Figure |2] shows that, after convolution 
with the LAT PSF, their PDFs are similar for a source that has a 0?5 radius typical of 
LAT-detected extended sources. To allow a valid comparison between the Gaussian and the 
uniform disk models, we define the source size as the radius containing 68% of the intensity 
(r 68 ). By direct integration, we find 

r 68, Gaussian =1-51(7, (9) 

r 6 8,disk =0.82a, (10) 

169 where a is defined in Equation [7] and Equation [8] respectively. For the example above, 

i/o r 68 = 0?5 so o- disk = 0.61° and cr Gaussian = 0.33°. 

171 For sources that are comparable in size to the PSF, the differences in the PDF for 

172 different spatial models are lost in the noise and the LAT is not sensitive to the detailed 

173 spatial structure of these sources. In section 13. 3[ we perform a dedicated Monte Carlo 

174 simulation that shows there is little bias due to incorrectly modeling the spatial structure 

175 of an extended source. Therefore, in our search for extended sources we use only a radially- 

176 symmetric uniform disk spatial model. Unless otherwise noted, we quote the radius to the 

177 edge (a) as the size of the source. 



3. Validation of the TS Distribution 

3.1. Point-like Source Simulations Over a Uniform Background 

We tested the theoretical distribution for TS cxt to evaluate the false detection proba- 
bility for measuring sour c e exte nsion. To do so, we tested simulated point-like sources for 
extension. iMattox et al.l (jl996[ ) discuss that the TS distribution for a likelihood-ratio test 



on the existence of a source at a given position is 

P(TS) = |(x?(TS) + 5(TS)), (11) 

where P(TS) is the probability density to get a particular value of TS, xl is the chi-squared 
distribution with one degree of freedom, and S is the Dirac delta function. The particular 
form of Equation [11] is due to the null hypothesis (source flux $ = 0) residing on the edge 
of parameter space and the model hypothesis adding a single degree of freedom (the source 
flux). It leads to the often quoted result y/TS = a, where a here refers to the significance of 
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185 the detection. It is plausible to expect a similar distribution for the TS in the test for source 

186 extension since the same conditions apply (with the source flux $ replaced by the source 

187 radius r and r < being unphysical). To verify Equation [TT], we evaluated the empirical 

188 distribution function of TS cxt computed from simulated sources. 

189 We simulated point-like sources with various spectral forms using the LAT on-orbit 

190 simulation tool gtobssinjfl and fit the sources with pointlike using both point-like and ex- 

191 tended source hypotheses. These point-like sources were simulated with a power-law spectral 

192 model with integrated fluxes above 100 MeV ranging from 3 x 10 -9 to 1 x 10~ 5 phcm~ 2 s" 1 

193 and spectral indices ranging from 1.5 to 3. These values were picked to represent typical 

194 parameters of LAT-detected sources. The point-like sources were simulated on top of an 

195 isotropic background with a power-law spectral model with integrated flux above 100 MeV 
we of 1.5 x 10~ 5 phcm -2 s _1 sr _1 and spectra l index 2.1. This was t aken to be the same as the 



197 isotropic spectrum measured by EGRET (ISreekumar et al.lll998l ). This spectrum is compa 



198 rable to the high-latitude background intensity seen by the LAT. The Monte Carlo simulation 

199 was performed over a one- year observation period using a representative spacecraft orbit and 

200 liveti me. The reconstruct ion was performed using the P7_V6 Source class event selection and 



201 IRFs (lAbdo et al.ll2012al ). For each significantly detected point-like source (TS > 25), we 

202 used point like to fit the source as an extended source and calculate TS ext . This entire 

203 procedure was performed twice, once fitting in the 1 GeV to 100 GeV energy range and once 

204 fitting in the 10 GeV to 100 GeV energy range. 

205 For each set of spectral parameters, ~ 20, 000 statistically independent simulations 

206 were performed. For lower-flux spectral models, many of the simulations left the source 

207 insignificant (TS < 25) and were discarded. Table [1] shows the different spectral models 

208 used in our study as well as the number of simulations and the average point-like source 

209 significance. The cumulative density of TS ext is plotted in Figures [3] and H] and compared to 

210 the Xi/2 distribution of Equation [TT1 

2n Our study shows broad agreement between simulations and Equation [TTJ To the extent 

212 that there is a discrepancy, the simulations tended to produce smaller than expected values 

213 of TS cxt which would make the formal significance conservative. Considering the distribution 

214 in Figures [3] and HI the choice of a threshold TS ex t set to 16 (corresponding to a formal 4a 

215 significance) is reasonable. 



3 gtobssim is distributed publicly by the FSSC. 
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216 3.2. Point-like Source Simulations Over a Structured Background 

217 We performed a second set of simulations to show that the theoretical distribution 
2ib for TScxt is still preserved when the point-like sources are present over a highly-structured 

219 diffuse background. Our simulation setup was the same as above except that the sources 

220 were simulated on top of and analyzed assuming the presence of the standard Galactic diffuse 

221 and isotropic background models used in 2FGL. In our simulations, we selected our sources 

222 to have random positions on the sky such that they were within 5° of the Galactic plane. 

223 This probes the brightest and most strongly contrasting areas of the Galactic background. 

224 To limit the number of tests, we selected only one flux level for each of the four spectral 

225 indices and we performed this test only in the 1 GeV to 100 GeV energy range. As described 

226 below, the fluxes were selected so that TS ~ 50. We do not expect to be able to spatially 

227 resolve sources at lower fluxes than these, and the results for much brighter sources are less 
22B likely to be affected by the structured background. 

Because the Galactic diffuse emission is highly structured with strong gradients, the 
point-source sensitivity can vary significantly across the Galactic plane. To account for this, 
we scaled the flux (for a given spectral index) so that the source always has approximately 
the same signal-to-noise ratio: 

F "> fiGc > x (H)"' < i2 > 

229 Here, x is the position of the simulated source, F is the integral flux of the source from 100 

230 MeV to 100 GeV, F(GC) is the same quantity if the source was at the Galactic center, B 

231 is the integral of the Galactic diffuse and isotropic emission from 1 GeV to 100 GeV at the 

232 position of the source, and -B(GC) is the same quantity if the source was at the Galactic 

233 center. For the four spectral models, Table [T] lists F(GC) and the average value of TS. 

234 For each spectrum, we performed ~ 90, 000 simulations. Figure [5] shows the cumulative 

235 density of TS cxt for each spectrum. For small values of TS ext , there is good agreement 

236 between the simulations and theory. For the highest values of TS ext , there is possibly a small 

237 discrepancy, but the discrepancy is not statistically significant. Therefore, we are confident 

238 we can use TS cx t as a robust measure of statistical significance when testing LAT-detected 

239 sources for extension. 
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240 3.3. Extended Source Simulations Over a Structured Background 

241 We also performed a Monte Carlo study to show that incorrectly modeling the spatial 

242 extension of an extended source does not substantially bias the spectral fit of the source, 

243 although it does alter the value of the TS. To assess this, we simulated the spatially extended 

244 ring-type SNR W44. We selected W44 because it is the most significant exte nded source 



245 detect ed by the LAT that has a non-radially symmetric photon distribution ( lAbdo et al. 



2010gj ). 



W44 was simulated with a power-law spectral model with an integral flux of 7.12 x 10 -8 
phcm _2 s~ 1 in the energy range from 1 GeV to 100 GeV and a spectral index of 2.66 (see 
Section [7]). 



W44 was simulated with the elliptical ring spatial model described in lAbdo et al.l ( I2010gh 
For reference, the ellipse has a semi-major axis of 0?3, a semi-minor axis of 0?19, a position 
angle of 147° measured East of celestial North, and the ring's inner radius is 75% of the 
outer radius. 



254 We used a simulation setup similar to that described in Section 13. 2\ but the simulations 

255 were over the 2-year interval of the 2FGL catalog. In the simulations, we did not include 

256 the finite energy resolution of the LAT to isolate any effects due to changing the assumed 

257 spatial model. The fitting code we use also ignores this energy dispersion and the potential 
25B bias introduced by this will be discussed in an upcoming paper by the LAT collaboration 



259 f lAbdo et al.ll2012al ). In total, we performed 985 independent simulations. 



260 The simulated sources were fit using a point-like spatial model, a radially-symmetric 

261 Gaussian spatial model, a uniform disk spatial model, an elliptical disk spatial model, and 

262 finally with an elliptical ring spatial model. We obtained the best fit spatial parameters 

263 using point like and, with these parameters, obtained the best fit spectral parameters using 

264 gtlike. 

265 Figure |6^ shows that the significance of W44 in the simulations is very large (TS ~ 3500) 

266 for a model with a point-like source hypothesis. Figure [6b shows that the significance of 

267 the spatial extension is also large (TS cx t ~ 250). On average TS ex t is somewhat larger when 

268 fitting the sources with more accurate spatial models. This shows that assuming an incorrect 

269 spatial model will cause the source's significance to be underestimated. Figure [6b shows that 

270 the sources were fit better when assuming an elliptical disk spatial model compared to a 

271 uniform disk spatial model (TS e iii P ticai disk — TSdisk ~ 30). Finally, Figure [6U shows that the 

272 sources were fit somewhat better assuming an elliptical ring spatial model compared to an 

273 elliptical disk spatial model (TS e iii P ticai ring — TS c iiipti ca i disk ~ 9). This shows that the LAT 

274 has some additional power to resolve substructure in bright extended sources. 
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275 Figure [T^i and Figure [7b clearly show that no significant bias is introduced by modeling 

276 the source as extended but with an inaccurate spatial model, while a point-like source model- 

277 ing results in a ~ 10% and ~ 0.125 bias in the fit flux and index, respectively. Furthermore, 
27B Figure [Tfc shows that the r 68 estimate of the extension size is very mildly biased (~ 10%) 

279 toward higher values when inaccurate spatial models are used, and thus represents a rea- 

280 sonable measurement of the true 68% containment radius for the source. For the elliptical 

281 spatial models, r 68 is computed by numeric integration. 

282 4. Extended Source Detection Threshold 

283 We calculated the LAT flux threshold to detect spatial extent. We define the detection 

284 threshold as the flux at which the value of TS cx t averaged over many statistical realizations 

285 is (TSext) — 16 (corresponding to a formal 4er significance) for a source of a given extension. 

286 We used a simulation setup similar to that described in Section 13. 1\ but instead of point- 

287 like sources we simulated extended sources with radially-symmetric uniform disk spatial 

288 models. Additionally, we simulated our sources over the two-year time range included in 

289 the 2FGL catalog. For each extension and spectral index, we selected a flux range which 

290 bracketed TS ext = 16 and performed an extension test for > 100 independent realizations of 

291 ten fluxes in the range. We calculated (TS ex t) = 16 by fitting a line to the flux and TS cxt 

292 values in the narrow range. 

293 Figure [8] shows the threshold for sources of four spectral indices from 1.5 to 3 and 

294 extensions varying from a = 0?1 to 2?0. The threshold is high for small extensions when 

295 the source is small compared to the size of the PSF. It drops quickly with increasing source 

296 size and reaches a minimum around 0?5. The threshold increases for large extended sources 

297 because the source becomes increasingly diluted by the background. Figure [8] shows the 

298 threshold using photons with energies between 100 MeV and 100 GeV and also using only 

299 photons with energies between 1 GeV and 100 GeV. Except for very large or very soft sources, 

300 the threshold is not substantially improved by including photons with energies between 100 

301 MeV and 1 GeV. This is also demonstrated in Figure [j] which shows TS ex t for the SNR 

302 IC 443 computed independently in twelve energy bins between 100 MeV and 100 GeV. For 

303 IC 443, which has a spectral index ~ 2.4 and an extension ~ 0?35, almost the entire increase 

304 in likelihood from optimizing the source extent in the model comes from energies above 1 

305 GeV. Furthermore, other systematic errors become increasingly large at low energy. For our 

306 extension search (Section [9]), we therefore used only photons with energies above 1 GeV. 

307 Figure M shows the flux threshold as a function of source extension for different back- 
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308 ground levels (lx, 10 x, and 100 x the nominal background), different spectral indices, and 

309 two different energy ranges (1 GeV to 100 GeV and 10 GeV to 100 GeV). The detection 

310 threshold is higher for sources in regions of higher background. When studying sources only 

311 at energies above 1 GeV, the LAT detection threshold (defined as the 1 GeV to 100 GeV flux 

312 at which (TS cxt ) = 16) depends less strongly on the spectral index of the source. The index 

313 dependence of the detection threshold is even weaker when considering only photons with 

314 energies above 10 GeV because the PSF changes little from 10 GeV to 100 GeV. Overlaid 

315 on Figure M are the LAT-detected extended sources that will be discussed in Sections [7] and 

316 [TUl The extension thresholds are tabulated in Table [2j 

317 Finally, Figure [10] shows the projected detection threshold of the LAT to extension 

318 with a 10 year exposure against 10 times the isotropic background measured by EGRET. 

319 This background is representative of the background near the Galactic plane. For small 

320 extended sources, the threshold improves by a factor larger than the square root of the 

321 relative exposures because the LAT is signal-limited at high energies where the present 

322 analysis is most sensitive. For large extended sources, the relevant background is over a 

323 larger spatial range and so the improvement is closer to a factor corresponding to the square 

324 root of the relative exposures that is caused by Poisson fluctuations in the background. 

325 5. Testing Against Source Confusion 

326 It is impossible to discriminate using only LAT data between a spatially extended source 

327 and multiple point-like sources separated by angular distances comparable to or smaller than 

328 the size of the LAT PSF. To assess the plausibility of source confusion for sources with 

329 TS cxt > 16, we developed an algorithm to test if a region contains two point-like sources. 

330 The algorithm works by simultaneously fitting in point like the positions and spectra of the 

331 two point-like sources. To help with convergence, it begins by dividing the source into two 

332 spatially coincident point-like sources and then fitting the sum and difference of the positions 

333 of the two sources without any limitations on the fit parameters. 

After simultaneously fitting the two positions and two spectra, we define TS2 P t s as twice 
the increase in the log of the likelihood fitting the region with two point-like sources compared 
to fitting the region with one point-like source: 



334 For the following analysis of LAT data, TS2 P t s was computed by fitting the spectra of the two 

335 point-like sources in gtlike using the best fit positions of the sources found by pointlike. 




(13) 



336 



TS2 P ts cannot be quantitatively compared to TS ext using a simple likelihood-ratio test to 
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337 evalu ate which model is significantly better because the models are not nested fjProtassov et al 



338 120021 ) . Even though the comparison of TS ext with TS2 P t s is not a calibrated test, TS cxt > 

339 TS2 P ts indicates that the likelihood for the extended source hypothesis is higher than for two 

340 point-like sources and we only consider a source to be extended if TS ext > TS 2pts . 



341 We considered using the Bayesian information criterion (BIC, ISchwarzl 119781 ) as an 

342 alternative Bayesian formulation for this test, but it is difficult to apply to LAT data because 

343 it contains a term including the number of data points. For studying 7-ray sources in LAT 

344 data, we analyze relatively large regions of the sky to better define the contributions from 

345 diffuse backgrounds and nearby point sources. This is important for accurately evaluating 

346 source locations and fluxes but the fraction of data directly relevant to the evaluation of the 

347 parameters for the source of interest is relatively small. 



As an alternative, we considered the Akaike information criterion test (AIC, lAkaike 



349 Il974l ). The AIC is defined as AIC = 2k — 2 log £, where k is the number of parameters in the 

350 model. In this formulation, the best hypothesis is considered to be the one that minimizes 

351 the AIC. The first term penalizes models with additional parameters. 

352 The two point-like sources hypothesis has three more parameters than the single ex- 

353 tended source hypothesis (two more spatial parameters and two more spectral parameters 

354 compared to one extension parameter), so the comparison AIC cx t < AIC2 P ts is formally 

355 equivalent to TS cxt + 6 > TS 2pts . Our criterion for accepting extension (TS cxt > TS 2p t s ) 

356 is thus equivalent to requesting that the AlC-based empirical support for the two point- 

357 like sources mo del be "considerably les s " than for the extended source model, following the 



358 classification by lBurnham &: Anderson! (120021 ) . 



359 We assessed the power of the TS ext > TS 2p t s test with a Monte Carlo study. We 

360 simulated one spatially extended source and fit it as both an extended source and as two 

361 point-like sources using pointlike. We then simulated two point-like sources and fit them 

362 with the same two hypotheses. By comparing the distribution of TS 2pts and TS cx t computed 

363 by pointlike for the two cases, we evaluated how effective the TS cxt > TS 2pts test is at 

364 rejecting cases of source confusion as well as how likely it is to incorrectly reject that an 

365 extended source is spatially extended. All sources were simulated using the same time range 

366 as in Section H] against a background 10 times the isotropic background measured by EGRET, 

367 representative of the background near the Galactic plane. 

368 We did this study first in the energy range from 1 GeV to 100 GeV by simulating 

369 extended sources of flux 4 x 10~ 9 ph cm -2 s~ x integrated from 1 GeV to 100 GeV and a power- 

370 law spectral model with spectral index 2. This spectrum was picked to be representative of 

371 the new extended sources that were discovered in the following analysis when looking in the 
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372 1 GeV to 100 GeV energy range (see Section [TO]) . We simulated these sources using uniform 

373 disk spatial models with extensions varying up to 1°. Figure ITTa shows the distribution of 

374 TS cxt and TS 2p t s and Figure [TTb shows the distribution of TS ex t — TS 2p t s as a function of the 

375 simulated extension of the source for 200 statistically independent simulations. 

376 Figure [T2k shows the same plot but when fitting two simulated point-like sources each 

377 with half of the flux of the spatially extended source and with the same spectral index as the 
37B extended source. Finally, Figure [T2"b shows the same plot with each point-like source having 

379 the same flux but different spectral indices. One point-like source had a spectral index of 

380 1.5 and the other an index of 2.5. These indices are representative of the range of indices of 

381 LAT-detected sources. 

382 The same four plots are shown in Figure [TTb . [TTH . [T2"b . and [T2H but this time when 

383 analyzing a source of flux 10~ 9 phcm -2 s _1 (integrated from 10 GeV to 100 GeV) only in the 

384 1 GeV to 100 GeV energy range. This flux is typical of the new extended sources discovered 

385 using only photons with energies between 10 GeV and 100 GeV (see Section [TO]) . 

386 Several interesting conclusions can be made from this study. As one would expect, 

387 TScxt — TS2 P ts is mostly positive when fitting the simulated extended sources. In the 1 GeV 

388 to 100 GeV analysis, only 11 of the 200 simulated extended sources had TS cx t > 16 but 

389 were incorrectly rejected due to TS2 P t s being greater than TS ext . In the 10 GeV to 100 GeV 

390 analysis, only 7 of the 200 sources were incorrectly rejected. From this, we conclude that 

391 this test is unlikely to incorrectly reject truly spatially extended sources. 

392 On the other hand, it is often the case that TS cx t > 16 when testing the two simulated 

393 point-like sources for extension. This is especially the case when the two sources had the 

394 same spectral index. Forty out of 200 sources in the 1 GeV to 100 GeV energy range and 

395 43 out of 200 sources in the 10 GeV to 100 GeV energy range had TS ext > 16. But in these 

396 cases, we always found the single extended source fit to be worse than the two point-like 

397 source fit. From this, we conclude that the TS cx t > TS2 P ts test is powerful at discarding 

398 cases in which the true emission comes from two point-like sources. 

399 The other interesting feature in Figure [TTk and [TTb is that for simulated extended 

400 sources with typical sizes (a ~ 0?5), one can often obtain almost as large an increase in 

401 likelihood fitting the source as two point-like sources (TS2 P t s ~ TS ext ). This is because 

402 although the two point-like sources represent an incorrect spatial model, the second source 

403 has four additional degrees of freedom (two spatial and two spectral parameters) and can 

404 therefore easily model much of the extended source and statistical fluctuations in the data. 

405 This effect is most pronounced when using photons with energies between 1 GeV and 100 

406 GeV where the PSF is broader. 
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w From this Monte Carlo study, we can see the limits of an analysis with LAT data of 

408 spatially extended sources. Section [37T1 showed that we have a statistical test that finds when 

409 a LAT source is not well described by the PSF. But this test does not uniquely prove that 

410 the emission originates from spatially extended emission instead of from multiple unresolved 

411 sources. Demanding that TS cx t > TS2 P t s is a powerful second test to avoid cases of simple 

412 confusion of two point-like sources. But it could always be the case that an extended source 

413 is actually the superposition of multiple point-like or extended sources that could be resolved 

414 with deeper observations of the region. There is nothing about this conclusion unique to 

415 analyzing LAT data, but the broad PSF of the LAT and the density of sources expected to 

416 be GeV emitters in the Galactic plane makes this issue more significant for analyses of LAT 

417 data. When possible, multiwavelength information should be used to help select the best 
4ib model of the sky. 



Test of 2LAC Sources 



420 For all following analyses of LAT data, we used the same two-year dataset that was 

421 used in the 2FGL catalog spanning from 2008 August 4 to 2010 August 1. We applied the 

422 same acceptance cu ts and we used the same P7_V6 Source class event selection and IRFs 

423 flAbdo et al.ll2012af ). When analyzing sources in pointlike, we used a circular 10° region 

424 of interest (ROI) centered on our source and eight energy bins per logarithmic decade in 

425 energy. When refitting the region in gtlike using the best fit spatial and spectral models 

426 from pointlike, we used the 'binned likelihood' mode of gtlike on a 14° x 14° ROI with 

427 a pixel size of 0?03. 

428 Unless explicitly mentioned, we used the same background model as 2FGL to represent 

429 the Galactic diffuse, isotropic, and Earth limb emission. To compensate for possible residuals 

430 in the diffuse emission model, the Galactic emission was scaled by a power-law and the 

431 normalization of the isotropic component was left free. Unless explicitly mentioned, we used 

432 all 2FGL sources within 15° of our source as our list of background sources and we refit the 

433 spectral parameters of all sources within 2° of the source. 

434 To validate our method, we tested LAT sources associated with AGN for extension. 

435 GeV emission from AGN is believed to originate from collimated jets. Therefore AGN are 

436 not expected to be spatially resolvable by the LAT and provide a good calibration source 

437 to demonstrate that our extension detection method does not misidentify point-like sources 

438 as being extended. We note that megapa rsec-scale 7-ray halos around AGNs have been 



439 hypothesized to be resolvable by the LAT (lAharonian et a 



440 has been discovered in the LAT data so far (INeronov et al 



20111 ) 



1994T) . However, no such halo 
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441 Following 2FGL, the LAT Collaboration published the Second LAT AGN Catalog 

442 (2LAC), a list of high latitude (\b\ > 10°) sources that had a high probability associa- 



te 



443 tion with AGN flAckermann et all 120111 ). 2LAC associated 1016 2FGL sources with AGN. 

444 To avoid systematic problems with AGN classification, we selected only the 885 AGN which 

445 made it into the clean AGN sub-sample defined in the 2LAC paper. An AGN association 

446 is considered clean only if it has a high probability of association P > 80%, if it is the only 

447 AGN associated with the 2FGL source, and if no analysis flags have been set for the source 
44B in the 2FGL catalog. These last two conditions are important for our analysis. Source con- 
449 fusion may look like a spatially extended source and flagged 2FGL sources may correlate 



with unmodeled structure in the diffuse emission. 



451 Of the 885 clean AGN, we selected the 733 of these 2FGL sources which were significantly 

452 detected above 1 GeV and fit each of them for extension. The cumulative density of TS cxt 

453 for these AGN is compared to the Xi/2 distribution of Equation [TT1 in Figure [TBI The TS cx t 

454 distribution for the AGN shows reasonable agreement with the theoretical distribution and 

455 no AGN was found to be significantly extended (TS cx t > 16). The observed discrepancy 

456 from the theoretical distribution is likely due to small systematics in our model of the LAT 

457 PSF and the Galactic diffuse emission (see Section [S]). The discrepancy could also in a few 

458 cases be due to confusion with a nearby undetected source. We note that the Monte Carlo 

459 study of section [37T1 effectively used perfect IRFs and a perfect model of the sky. The overall 

460 agreement with the expected distribution demonstrates that we can use TS ext as a measure 

461 of the statistical significance of the detection of the extension of a source. 

462 We note that the LAT PSF used in this study was determined empirically by fitting the 

463 distributions of gamma rays around bright AGN (see Section [8]). Finding that the AGN we 

464 test are not extended is not surprising. This validation analysis is not suitable to reject any 

465 hypotheses about the existence of megaparsec-scale halos around AGN. 
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7. Analysis of Extended Sources Identified in the 2FGL Catalog 

As further validation of our method for studying spatially extended sources, we re- 
analyzed the twelv e spatially extended sources which were included in the 2FGL catalog 



469 (iNolan et al.ll2012l ). Even though these sources had all been the subjects of dedicated analy- 

470 ses and separate publications, and had been fit with a variety of spatial models, it is valuable 

471 to show that these sources are significantly extended using our systematic method assuming 

472 radially-symmetric uniform disk spatial models. On the other hand, for some of these sources 

473 a uniform disk spatial model does not well describe the observed extended emission and so 

474 the dedicated publications by the LAT collaboration provide better models of these sources. 
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Six extended SNRs wer e included in the 2 FGL c atal og: W51C, IC 443 W28, W44 
the Cy gnus Loop, and W30 f kbdo et al.ll2009al . boiOhffl iKatagiri et aDboill ; kbdo et al. 
2012b|). Using photons with energies between 1 GeV and 100 GeV, our analysis significantly 
detected that these six SNRs are spatially extended. 

Two nearby satellite galaxies of the Milky Way the Large Magellanic Cloud (LMC) 
and the Small M agellanic Cloud (S MC) were included in the 2FGL catalog as spatially 
extended sources ( lAbdo et al.ll2010i| jbl). Their extensions were significantly detected using 
photons with energies between 1 GeV and 100 GeV. Our fit extensions are comparable to 
the published result, but we note that the previous LAT Collaboration publication on the 
LMC used a more complicated two 2D Gaussian surface brightness profile when fitting it 
( |Abdo et alihoiOih . 



Three PWNe, MSH 15-52, Vela X, and HESS J1825-137, were fit as extended sources 



in the 2FGL analysis ( lAbdo et al.l l2010al Jfl; iGrondin et al.ll201ll ). In the present analysis, 
HESS J1825— 137 was significantly detected using photons with energies between 10 GeV and 
100 GeV. To avoid confusion with the nearby bright pulsar PSR J1509— 5850, MSH 15— 52 
had to be analyzed at high energies. Using photons with energies above 10 GeV, we fit the 
extension of MSH 15—52 to be consistent with the published size but with TS cx t =6.6. 

Our analysis was unable to resolve Vela X which would have required first removing the 
pulsed photons from the Vela pulsar which was beyond the scope of this paper. Our analysis 
also failed to detect a significant extension for the Centaurus A Lobes becaus e the shape 
of the source is significantly different from a uniform radially-symmetric disk (lAbdo et al. 
2010ch . 



Our analysis of these sources is summarized in Table [3J This table includes the best 
fit positions and extensions of these sources when fitting them with a radially-symmetric 
uniform disk model. It also includes the best fit spectral parameters for each source. The 
positions and extensions of Vela X and the Centaurus A Lobes were taken from AbdoetaL 
and are included in this table for completeness. 



8. Systematic Errors on Extension 



We developed two criteria for estimating systematic errors on the extensions of the 
sources. First, we estimated a systematic error due to uncertainty in our knowledge of the 
LAT PSF. Before launch, the LAT PSF was determined by detector simulations which were 
verified in accelerator beam tests ( lAtwood et al.l l2009bl ). However, in-flight data revealed 
a discrepancy above 3 GeV in the PSF compared to the angular distribution of photons 
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from bright AGN (lAbdo et al.ll2012al ). Subsequently, the PSF was fit empirically to bright 
AGN and this empirical parameterization is used in the P7_V6 IRFs. To account for the 
uncertainty in our knowledge of the PSF, we refit our extended source candidates using the 
pre-flight Monte Carlo representation of the PSF and consider the difference in extension 
found using the two PSFs as a systema tic error on the extension of a source. The same 
approach was used in lAbdo et al.l ( 12010hl ). We believe that our parameterization of the PSF 
from bright AGN is substantially better than the Monte Carlo representation of the PSF so 
this systematic error is conservative. 

We estimated a second systematic error on the extension of a source due to uncertainty 
in our model of the Galactic diffuse emission by using an alternative approach to modeling 
the diffuse emission which takes as input templates calculated by GALPROP0 but then fits 
each template locally in the surrounding region. The particular GALPROP model that was 
used as input is de scribed in the analysis of the isotropic diffuse emission with LAT data 
( lAbdo et al.ll2010jl ). The intensities of various components of the Galactic diffuse emission 
were fitted individually using a spatial distribution predicted by the model. We considered 
separate contributions from cosmic-ray interactions with the molecular hydrogen, the atomic 
and ionized hydrogen, residual gas traced by dust ( jGrenier et al.l 120051 ). and the interstellar 
radiation field. We further split the contributions from interactions with molecular and 
atomic hydrogen to the Galactic diffuse emission according to the distance from the Galactic 
center in which they are produced. Hence, we replaced the standard diffuse emission model 
by 18 individually fitted templates to describe individual components of the diffuse emission. 
A similar crossche ck was used in an analysis of RX J1713.7— 3946 by the LAT Collaboration 
f lAbdo et al. l boilh . 



It is not expected that this diffuse model is superior to the standard LAT model obtained 
through an all-sky fit. However, adding degrees of freedom to the background model can 
remove likely spurious sources that correlate with features in the Galactic diffuse emission. 
Therefore, this tests systematics that may be due to imperfect modeling of the diffuse emis- 
sion in the region. Nevertheless, this alternative approach to modeling the diffuse emission 
does not test all systematics related to the diffuse emission model. In particular, because 
the alternative approach uses the same underlying gas maps, it is unable to be used to assess 
systematics due to insufficient resolution of the underlying maps. Structure in the diffuse 
emission that is not correlated with these maps will also not be assessed by this test. 



4 GALPROP is a software package for calculating the Galactic 7-ray emission based on a model of cosmic- 
ray propagation in the Galaxy and maps of the dist ributions of the components of the interstellar medium 
(jStrong fc Moskalenkolll998t IVladimirov et al.ll2011l ). See also http://galprop.stanford.edu/ for details. 
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540 We do not expect the systematic error due to uncertainties in the PSF to be correlated 

541 with the systematic error due to uncertainty in the Galactic diffuse emission. Therefore, the 

542 total systematic error on the extension of a source was obtained by adding the two errors in 

543 quadrature. 

544 There is another systematic error on the size of a source due to issues modeling nearby 

545 sources in crowded regions of the sky. It is beyond the scope of this paper to address this 

546 systematic error. Therefore, for sources in crowded regions the systematic errors quoted in 

547 this paper may not represent the full set of systematic errors associated with this analysis. 

54B 9. Extended Source Search Method 

549 Having demonstrated that we understand the statistical issues associated with analyzing 

550 spatially extended sources (Section 13.11 and E]) and that our method can correctly analyze 

551 the extended sources included in 2FGL (Section [7]), we applied this method to search for 

552 new spatially extended GeV sources. The data and general analysis setting is as described 

553 in Section [61 

554 Ideally, we would apply a completely blind and uniform search that tests the extension 

555 of each 2FGL source in the presence of all other 2FGL sources to find a complete list of 

556 all spatially extended sources. As our test of AGN in Section [6] showed, at high Galactic 

557 latitude where the source density is not as large and the diffuse emission is less structured, 

558 this method works well. 

559 But this is infeasible in the Galactic plane where we are most likely to discover new 

560 spatially extended sources. In the Galactic plane, this analysis is challenged by our imperfect 

561 model of the diffuse emission and by an imperfect model of nearby sources. The Monte 

562 Carlo study in Section [5] showed that the overall likelihood would greatly increase by fitting 

563 a spatially extended source as two point-like sources so we expect that spatially extended 

564 sources would be modeled in the 2FGL catalog as multiple point-like sources. Furthermore, 

565 the positions of other nearby sources in the region close to an extended source could be 
see biased by not correctly modeling the extension of the source. The 2FGL catalog contains a 

567 list of sources significant at energies above 100 MeV whereas we are most sensitive to spatial 

568 extension at higher energies. We therefore expect that at higher energies our analysis would 

569 be complicated by 2FGL sources no longer significant and by 2FGL sources whose positions 

570 were biased by diffuse emission at lower energies. 

571 To account for these issues, we first produced a large list of possibly extended sources 

572 employing very liberal search criteria and then refined the analysis of the promising can- 
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573 didates on a case by case basis. Our strategy was to test all point-like 2FGL sources for 

574 extension assuming they had a uniform radially-symmetric disk spatial model and a power- 

575 law spectral model. Although not all extended sources are expected to have a shape very 

576 similar to a uniform disk, Section 12.41 showed that for many spatially extended sources the 

577 wide PSF of the LAT and limited statistics makes this a reasonable approximation. On the 

578 other hand, choosing this spatial model biases us against finding extended sources that are 

579 not well described by a uniform disk model such as shell-type SNRs. 

580 Before testing for extension, we automatically removed from the background model all 

581 other 2FGL sources within 0?5 of the source. This distance is somewhat arbitrary, but was 

582 picked in hopes of finding extended sources with sizes on the order of ~ 1° or smaller. On the 

583 other hand, by removing these nearby background sources we expect to also incorrectly add 

584 to our list of extended source candidates point-like sources that are confused with nearby 

585 sources. To screen out obvious cases of source confusion, we performed the dual localization 

586 procedure described in Section [5]to compare the extended source hypothesis to the hypothesis 

587 of two independent point-like sources. 

588 As was shown in Section HI little sensitivity is gained by using photons with energies 

589 below 1 GeV. In addition, the broad PSF at low energy makes the analysis more susceptible 

590 to systematic errors arising from source confusion due to nearby soft point-like sources and 

591 by uncertainties in our modeling of the Galactic diffuse emission. For these reasons, we 

592 performed our search using only photons with energies between 1 GeV and 100 GeV. 

593 We also performed a second search for extended sources using only photons with en- 

594 ergies between 10 GeV and 100 GeV. Although this approach tests the same sources, it is 

595 complementary because the Galactic diffuse emission is even less dominant above 10 GeV 

596 and because source confusion is less of an issue. A similar procedure wa s used to detect the 



597 spatial extensions o f MSH 15-52 and HESS J1825-137 with the LAT flAbdo et al.ll2010a 



Grondin et all 1201 If ) 



599 When we applied this test to the 1861 point-like sources in the 2FGL catalog, our 

eoo search found 117 extended source candidates in the 1 GeV to 100 GeV energy range and 11 

eol extended source candidates in the 10 GeV to 100 GeV energy range. Most of the extended 

602 sources found above 10 GeV were also found above 1 GeV and in many cases multiple nearby 

603 point-like sources were found to be extended even though they fit the same emission region. 

604 For example, the sources 2FGL J1630.2— 4752, 2FGL J1632.4— 4753c 2FGL J1634.4— 4743c, 

605 and 2FGL J1636.3— 4740c were all found to be spatially extended in the 10 GeV to 100 GeV 

606 energy range even though they all fit to similar positions and sizes. For these situations, we 

607 manually discarded all but one of the 2FGL sources. 
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eos Similarly, many of these sources were confused with nearby point-like sources or influ- 

609 enced by large-scale residuals in the diffuse emission. To help determine which of these fits 

610 found truly extended sources and when the extension was influenced by source confusion and 
en diffuse emission, we generated a series of diagnostic plots. For each candidate, we generated 

612 a map of the residual TS by adding a new point-like source of spectral index 2 into the 

613 region at each position and finding the increase in likelihood when fitting its flux. Figure 

614 shows this map around the most significantly extended source IC 443 when it is modeled 

615 both as a point-like source and as an extended source. The residual TS map indicates that 

616 the spatially extended model for IC 443 is a significantly better description of the observed 

617 photons and that there is no TS > 25 residual in the region after modeling the source as 
6ib being spatially extended. We also generated plots of the sum of all counts within a given 

619 distance of the source and compared them to the model predictions assuming the emission 

620 originated from a point-like source. An example radial integral plot is shown for the ex- 

621 tended source IC 443 in Figure [TJ For each source, we also made diffuse-emission-subtracted 

622 smoothed counts maps (shown for IC 443 in Figured]). 

623 We found by visual inspection that in many cases our results were strongly influenced by 

624 large-scale residuals in the diffuse emission and hence the extension measure was unreliable. 

625 This was especially true in our analysis of sources in the 1 GeV to 100 GeV energy range. 

626 An example of such a case is 2FGL J1856.2+0450c analyzed in the 1 GeV to 100 GeV energy 

627 range. Figure [15] shows a diffuse-emission-subtracted smoothed counts map for this source 
62B with the best fit extension of the source overlaid. There appear to be large-scale residuals in 

629 the diffuse emission in this region along the Galactic plane. As a result, 2FGL J1856.2+0450c 

630 is fit to an extension of ~ 2° and the result is statistically significant with TS cx t =45.4. 

631 However, by looking at the residuals it is clear that this complicated region is not well 

632 modeled. We manually discard sources like this. 

633 We only selected extended source candidates in regions that did not appear dominated 

634 by these issues and where there was a multiwavelength counterpart. Because of these sys- 

635 tematic issues, this search can not be expected to be complete and it is likely that there are 

636 other spatially extended sources that this method missed. 

637 For each candidate that was not biased by neighboring point-like sources or by large-scale 

638 residuals in the diffuse emission model, we improved the model of the region by deciding on a 

639 case by case basis which background point-like sources should be kept. We kept in our model 

640 the sources that we believed represented physically distinct sources and we removed sources 

641 that we believed were included in the 2FGL catalog to compensate for residuals induced by 

642 not modeling the extension of the source. Soft nearby point-like 2FGL sources that were 

643 not significant at higher energies were frozen to the spectras predicted by 2FGL. When 
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644 deciding which background sources to keep and which to remove, we used multiwavelength 

645 information about possibly extended source counterparts to help guide our choice. For each 

646 extended source presented in Section [TD], we describe any modifications from 2FGL of the 

647 background model that were performed. In Table El we summarize the sources in the 2FGL 

648 catalog that we have concluded here correspond to residuals induced by not modeling the 

649 extensions of nearby extended sources. 

650 The best fit positions of nearby point-like sources can be influenced by the extended 

651 source and vice versa. Similarly, the best fit positions of nearby point-like sources in the 2FGL 

652 catalog can be biased by systematic issues at lower energies. Therefore, after selecting the list 

653 of background sources, we iteratively refit the positions and spectra of nearby background 

654 sources as well as the positions and extensions of the analyzed spatially extended sources until 

655 the overall fit converged globally. For each extended source, we will describe the positions 

656 of any relocalized background sources. 

657 After obtaining the overall best fit positions and extensions of all of the sources in 

658 the region using point like, we refit the spectral parameters of the region using gtlike. 

659 With gtlike, we obtained a second measure of TS ex t- We only consider a source to be 

660 extended when both pointlike and gtlike agree that TS cx t > 16. We further required 
eel that TS cxt > 16 using the alternative approach to modeling the diffuse emission presented 

662 in Section [81 We then replaced the spatially extended source with two point-like sources 

663 and refit the positions and spectra of the two point-like sources to calculate TS 2pts . We only 

664 consider a source to be spatially extended, instead of being the result of confusion of two 
ees point-like sources, if TS cx t > TS2 P ts- As was shown in Section [51 this test is fairly powerful 
eee at removing situations in which the emission actually originates from two distinct point- 
667 like sources instead of one spatially extended source. On the other hand, it is still possible 
ees that longer observations could resolve additional structure or new sources that the analysis 

669 cannot currently detect. Considering the very complicated morphologies of extended sources 

670 observed at other wavelengths and the high density of possible sources that are expected 

671 to emit at GeV energies, it is likely that in some of these regions further observations will 

672 reveal that the emission is significantly more complicated than the simple radially-symmetric 

673 uniform disk model that we assume. 



10. New Extended Sources 



675 Nine extended sources not included in the 2FGL catalog were found by our extended 

676 source search. Two of these hav e been previously stu died in dedicated publications: RX J1713.7— 3946 

677 and Vela Jr. fjAbdo et al.ll201ll ; iTanaka et al.ll201ll ). Two of these sources were found when 



-24- 



67B using photons with energies between 1 GeV and 100 GeV and seven were found when using 

679 photons with energies between 10 GeV and 100 GeV. For the sources found at energies above 

680 10 GeV, we restrict our analysis to higher energies because of the issues of source confusion 

681 and diffuse emission modeling described in Section The spectral and spatial properties 

682 of these nine sources are summarized in Table H] and the results of our investigation of sys- 

683 tematic errors are presented in Table |5j Table [5] also compares the likelihood assuming the 

684 source is spatially extended to the likelihood assuming that the emission originates from 
ess two independent point-like sources. For these new extended sources, TS ext > TS 2p t s so we 
ese conclude that the GeV emission does not originate from two physically distinct point-like 
687 sources (see Section [5]). Table [5] also includes the results of the extension fits using variations 
ess of the PSF and the Galactic diffuse model described in Section [HJ There is good agreement 
689 between TS ext and the fit size using the standard analysis, the alternative approach to mod- 
ego eling the diffuse emission, and the alternative PSF. This suggests that the sources are robust 
69i against mis-modeled features in the diffuse emission model and uncertainties in the PSF. 



10.1. 2FGL J0823. 0—4246 



693 
694 
695 
696 
697 
698 
699 
700 
701 

702 
703 
704 
705 
706 
707 
708 
709 



2FGL J0823. 0—4246 was found by our search to be an extended source candidate in the 
1 GeV to 100 GeV energy range and is spatially coincident with the SNR Puppis A. Figure ITB1 
shows a counts map of this source. There are two nearby 2FGL sources 2FGL J0823.4— 4305 
and 2FGL J0821. 0—4254 that are also coincident with the SNR but that do not appear 
to represent physically distinct sources. We conclude that these nearby point-like sources 
were included in the 2FGL catalog to compensate for residuals induced by not modeling the 
extension of this source and removed them from our model of the sky. After removing these 
sources, 2FGL J0823.0-4246 was found to have an extension a = 0?37 ± 0?03 stat ± 0?02 sys 
with TS cxt = 48.0. Figure [T7l shows the spectrum of this source. 

Puppis A has been st udied in detail in radio ( jCastelletti et al. hood ), and X-ray ( IPetre et al 



19961 : iHwang et all 120081 ). The fit extension of 2FGL J0823. 0—4246 ma tches well the size of 



Pupp is A in X-ray. The distance of Puppis A was estimated at 2.2 kpc ( iReynoso et al.lll995 



20031 ) and leads to a 1 GeV to 100 GeV luminosity of ~ 3 x 10 34 ergs s \ No molecular 



clouds have been observed directl y adjacent to Puppis A (IParon et al.l 120081 ). similar to the 
LAT-detected Cygnus Loop SNR ( jKatagiri et al.l 1201 ll ). The luminosity of Pup pis A is also 



smaller than tha t of other SNRs believed to interact with molecular clouds ( lAbdo et al. 
2009ai l2010hl l3B lAbdo et al.llioioh . 
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10.2. 2FGLJ0851. 7-4635 



711 2FGL J085 1.7— 4635 was found by our search to be an extended source candidate in the 

712 10 GeV to 100 GeV energy range and is spatially coincid ent with the SNR V ela Jr. This 

713 source was recently studied by the LAT Collaboration in iTanaka et al.l ( 1201 ll ). Figure [T8l 

714 shows a counts map of the source. Overlaid on Figure [TS] are TeV contours of Vela Jr. 

715 ( jAharonian et al.l l2007al ). There are three point-like 2FGL sources 2FGL J0848.5— 4535, 

716 2FGL J0853.5— 4711, and 2FGL J0855.4— 4625 which correlate with the multiwavelength 

717 emission of this SNR but do not appear to be physically distinct sources. They were most 
7ib likely included in the 2FGL catalog to compensate for residuals induced by not modeling 

719 the extension of Vela Jr. and were removed from our model of the sky. 

720 With this model of the background, 2FGL J0851.7— 4635 was found to have an extension 



of a = 1?15 ± 0?08 stat ± 0?02 sys with TS 



Jcxt 



86.8. The LAT size matches well the TeV 

722 morphology of Vela Jr. While fitting the extension of 2FGL J0851.7— 4635, we iteratively 

723 relocalized the position of the nearby point-like 2FGL source 2FGL J0854.7— 4501 to (/, b) = 

724 (266?24, 0?49) to better fit its position at high energies. 



725 10.3. 2FGL J1615. 0—5051 

726 2FGL J1615. 0—5051 and 2FGL J1615. 2-5138 were both found to be extended source 

727 candidates in the 10 GeV to 100 GeV energy range. Because they are less than 1° away 
72B from each other, they needed to be analyzed simultaneously. 2FGL J1615. 0—5051 is spa- 

729 tially coincident with the extended TeV source HESS J1616— 508 and 2FGL J1615.2— 5138 is 

730 coincident with the extended TeV source HESS J1614— 518. Figure [T9l shows a counts map 

731 of these sources and ov erlays the TeV contours of HESS J1616— 508 and HESS J1614— 518 



732 f lAharonian et al.l 120061 ). The figure shows that the 2FGL source 2FGL J1614.9— 5212 is 

733 very close to 2FGL J1615.2— 5138 and correlates with the same extended TeV source as 

734 2FGL J1615.2— 5138. We concluded that this source was included in the 2FGL catalog to 

735 compensate for residuals induced by not modeling the extension of 2FGL J1615.2— 5138 and 

736 removed it from our model of the sky. 

737 With this model of the sky, we iteratively fit the extensions of 2FGL J1615. 0—5051 
73b and 2FGL J1615.2— 5138. 2FGL J1615. 0—5051 was found to have an extension a = 0?32 ± 

739 0?04 stat ± 0?01 sys and TS cxt =16.7. 

740 The TeV counterpart of 2FGL J1615. 0—5051 wa s fit with a radially-sym metric Gaussian 



74i surface brightness profile with a = 0?136 ± 0?008 ( jAharonian et al.l 120061 ). This TeV size 



742 corresponds to a 68% containment radius of tq§ = 0?21 ± 0?01, comparable to the LAT size 
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743 r 68 = 0?26±0?03. Figure [20] shows that the spectrum of 2FGL J1615.0-5051 at GeV energies 

744 connects to the spectrum of HESS J1616— 508 at TeV energies. 



745 
746 
747 
748 
749 
750 
751 
752 
753 
754 
755 
756 
757 
758 



HESS J 16 16-508 is located in the region of two SNRs RCW103 (G332.4-04) and Kes 3 2 
(G332. 4+0.1) but is not spatially coincident with either of them ( lAharonian et al.ll2006l ). 
HESS J1616-508 is near three pulsars PSR J1614-5048, PSR J1616-5109, and PSR J1617-5055. 

Only PSR J161 7— 5055 is energetically capable of 



( iTorii et al.l Il998l ; lLandi et al 



2007a) 



powering the TeV emission and lAharonian et al.l (120061 ) speculated that HESS J1616— 508 
could be a PWN powered by this young pulsar. Because HESS J1616— 508 is 9' away from 
PSR J1617— 5055, this would require an asymmetric X-ray PWNe to power the TeV emis- 
sion. Chandra ACIS observations revealed an underluminous PWN of size ~ 1' around the 
pulsar that was not ori ented towards the TeV emission, rendering this association uncertain 
( IKargaltsev et al.ll2009l ) . No other promising co unterparts were observed at X-ray and soft 7- 
ray energies by Suzaku (jMatsumoto et al.ll2007l). Swift/XKT, IBIS/ISGRBI, BeppoSAX and 



XMM-Newton flLandi et al.l l2007al ). 



Kargaltsev et al.l (120091 ) discovered additional diffuse 



emission towards the center of HESS J1616— 508 using archival radio and infared observa- 
tions. Deeper observations will likely be necessary to understand this 7-ray source. 



759 10.4. 2FGLJ1615. 2-5138 

760 2FGL J1615.2-5138 was found to have an extension a = 0?42 ± 0?04 stat ± 0.02 sys with 

761 TScxt = 46.5. To test for the possibility that 2FGL J1615.2— 5138 is not spatially extended 

762 but instead composed of two point-like sources (one of them represented in the 2FGL catalog 

763 by 2FGL J1614.9— 5212), we refit 2FGL J1615.2— 5138 as two point-like sources. Because 

764 TS 2p ts = 35.1 is less than TS cxt = 46.5, we conclude that this emission does not originate 

765 from two closely-spaced point-like sources. 

766 2FGL J1615.2— 5138 is spatially coincident with the extended TeV source HESS J1614— 51< 

767 H.E.S.S. measured a 2D Gaussian extension of a = 0?23 ± 0?02 and a = 0?15 ± 0?02 

768 in the semi-major and semi-minor axis. This corresponds to a 68% containment size of 

769 r 68 = 0?35 ± 0?03 and 0?23 ± 0?03, consistent with the LAT size r 68 = 0?34 ± 0?03. Figure M 

770 shows that the spectrum of 2FGL J1615.2— 5138 at GeV energies connects to the spectrum 

771 of HESS J1614— 518 at TeV energies. Further data collected by H.E.S.S. in 2007 resolve a 

772 double peaked structure at TeV energies but no spectral variation across this source, suggest - 



773 ing that the emission is not the confusion of physically separate sources (IRowell et al.ll2008l ). 



This double peaked structure is also hinted at in the LAT counts map in Fig ure [191 but is not 



775 very significant. The TeV source was also detected by CANGAROO-III ( IMizukami et al. 

776 bonh . 
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There are five nearby pulsars, but none are luminous eno ugh to provide the energy 

HESSJ1614-518 is spa - 



output required to power the 7-ray emission f lRowel 



et al 



2008 



tially coincident with a young open cluster Pismis 22 (lLandi et al.ll2007bl ; lRowell et al.ll2008l ). 
Suzaku detected two promising X-ray candidates. Source A is an extended source consistent 
with the peak of HESS J 16 14— 518 and source B coincident with Pismis 22 and tow ards the 
center but in a relatively dim region of HESS J1614— 518 (IMatsumoto et al.l 120081 ). Three 
hypotheses have been presented to explain this emission: either source A is an SNR powering 
the 7-ray emission; source A is a PWN powered by an undiscovered pulsar in either source A 
or B; and fin ally that the emission may arise from hadronic acceleration in the stellar winds 
of Pismis 22 jMizukami et alJboilh . 



10.5. 2FGL J1627.0— 2425c 



0?42±0?05 sta ,± 



2FGL J1627.0— 2425c was found by our search to have an extension o 
6 sys with TS cxt = 32.4 using photons with energies between 1 GeV and 100 GeV. Figure 
shows a counts map of this source. 

791 This source is in a region of remarkably complicated diffuse emission. Even though it is 

792 16° from the Galactic plane, this source is on top of the core of the Ophiuchus molecular cloud 

793 which contains massive star-forming regions that are bright in infrared. The region also has 

794 abundant molecular and atomic gas traced by CO and H I and significant dark gas found 

795 only by its association with dust emission (IGrenier et al.l 120051 ). Embedded star- forming 

796 regions make it even more challenging to measure the column density of dust. Infared and 

797 CO (7=1—7-0) c ontours are overlaid on Figure |2"T1 and sh ow good spatial correlation with 

798 the GeV emission (jYoung et al.lll986l ; Ide Geus et al.lll990l ). This source might represent 7- 

799 ray emission from the interactions of cosmic rays with interstellar gas which has not been 

800 accounted for in the LAT diffuse emission model. 



10.6. 2FGL J1632.4— 4753c 



802 2FGL J1632.4— 4753c was found by our search to be an extended source candidate in 

803 the 10 GeV to 100 GeV energy range but is in a crowded region of the sky. It is spa- 

804 tially coincident with the TeV source HESS J1632— 478. Fig ure [22b shows a counts map of 

805 this source and overlays TeV contours of HESS J1632— 478 (lAharonian et al.l 120061 ). There 

806 are six nearby point-like 2FGL sources that appear to represent physically distinct sources 

807 and were included in our background model: 2FGL J1630.2— 4752, 2FGL J1631.7— 4720c, 
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808 2FGL J1632.4-4820c, 2FGL J1635.4-4717c, 2FGL J1636.3-4740c, and 2FGL J1638.0-4703c. 

809 On the other hand, one point-like 2FGL source 2FGL J1634.4— 4743c correlates with the ex- 
sio tended TeV source and at GeV energies does not appear physically separate. It is very close 
8n to the position of 2FGL J1632.4— 4753c and does not show spatially separated emission in the 

812 observed photon distribution. We therefore removed this source from our model of the back- 

813 ground. Figure [22b shows the same region with the background sources subtracted. With 

814 this model, 2FGL J1632.4-4753c was found to have an extension a = 0?35±0?04 stat ±0?02 sys 
sis with TS ext = 26.9. While fitting the extension of 2FGL J1632.4— 4753c, we iteratively relo- 

816 calized 2FGL J1635.4— 4717c to (l,b) = (337?23, 0?35) and 2FGL J1636.3— 4740c to (l,b) = 

817 (336?97,-0?07). 

sis H.E.S.S measured an extension of a = 0?21 ± 0?05 and 0?06 ± 0?04 along the semi- 

819 major and semi-minor axes when fitting HESS J1632— 478 with an elliptical 2D Gaussian 

820 surface brightness profile. This corresponds to a 68% containment size r68 = 0?31 ± 0?08 

821 and 0?09 ± 0?06 along the semi-major and semi-minor axis, consistent with the LAT size 

822 r 68 = 0?29 ± 0?04. Figure M shows that the spectrum of 2FGL J1632.4— 4753c at GeV 

823 energies connects to the spectrum of HESS J1632— 478 at TeV energies. 



Aharonian et al.l ( 120061 ) argued that HESS J1632— 478 is positionally coincident with the 
hard X-ray source IGR J1632-4751 observed by ASCA, INTEG RAL, and XMM-Newton 
( ISugizaki et al.l l2001t iTomsick et al.l l2003t iRodriguez et al.l 120031 ) , but this source is sus- 
pected to be a Galactic X-Ray Binary so the 7-ray extension disfavors the association. Fur- 
ther observations by XMM-Newton discovered point-like emission coincident w ith the peak 
of th e H.E.S.S. source surrounded by extended emission of size ~ 32" x 15" f lBalbo et al. 



20101). They found in archival MGPS-2 data a spatially coincident extended radio source 



( jMurphy et al.ll2007l ) and argued for a single synchrotron and inverse Compton process pro- 
ducing the radio, X-ray, and TeV emission, likely due to a PWN. The increased size at TeV 
energies compared to X-ray energies has previously been observed in se veral aging PWNe 



i ncluding HESS J1825-137 riGaem 



(lAharonian et al.ll2006l ; iFunk et al. 



er et al.ll2003l : lAharonian et al.ll2006bh. HESS J1 640-465 

2007). and Vela X ( IMarkwardt fc Ogelmanlll995l ; lAharonian et al 



2006cJ) and can be explained by different synchrotron cooling times for the electrons that 



produce X-rays and 7-rays. 



10.7. 2FGL J1712.4— 3941 



839 2FGL J1712.4— 3941 was found by our search to be spatially extended using photons 

840 with energies between 1 GeV and 100 GeV. This source is spatially coincide nt with the 

841 SNR RX J1713.7— 3946 and was recently studied by the LAT Collaboration in lAbdo et al. 
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( 120111 ). To avoid issues related to uncertainties in the nearby Galactic diffuse emission at 
lower energy, we restricted our analysis only to energies above 10 GeV. Figure [23] shows a 
smoothed counts map of the source. Ab ove 10 GeV, the GeV em ission nicely correlates with 
the TeV contours of RX J1713.7-3946 (Uharonian et alil2007bh and 2FGL J1712.4-3941 fit 
to an extension o = 0?56 ± 0?04 stat ± 0?02 sys with TS mt = 38.5. 



847 10.8. 2FGL J1837.3— 0700c 

848 2FGL J1837.3— 0700c was found by our search to be an extended source candidate in 

849 the 10 GeV to 100 GeV energy range and is spatially coincident with the TeV source 

850 HESS J1837— 069. This source is in a complicated region. Figure |24"r sh ows a smoothed 



85i count s map of the region and overlays the TeV contours of HESS J 1837— 069 (lAharonian et al. 



852 |2006j). There are two very nearby point-like 2FGL sources, 2FGL J1836.8-0623c and 2FGL J1839.3-0558c, 

853 that clearly represent distinct sources. On the other hand, there is another source 2FGL J1835.5— 0649 

854 located between the three sources that appears to correlate with the TeV morphology of 

855 HESS J1837— 069 but at GeV energies does not appear to represent a physically distinct 

856 source. We concluded that this source was included in the 2FGL catalog to compensate 

857 for residuals induced by not modeling the extension of this source and removed it from our 

858 model of the sky. Figure I2~4b shows a counts map of this region after subtracting these 

859 background sources. After removing 2FGL J1835.5— 0649, we tested for source confusion 

860 by fitting 2FGL J1837.3— 0700c instead as two point-like sources. Because TS2 P ts = 10.8 is 

861 less than TS cx t = 18.5, we conclude that this emission does not originate from two nearby 

862 point- like sources. 

863 With this model, 2FGL J1837.3— 0700c was found to have an extension a = 0?33 ± 

864 0?07 sta t ± 0?05 sys . While fitting the extension of 2FGL J1837.3— 0700c, we iteratively relo- 
ses calized the two closest background sources along with the extension of 2FGL J1837.3— 0700c 
see but their positions did not significantly change. 2FGL J1834.7— 0705c moved to (7,5) = 
867 (24?77, 0?50), 2FGL J1836.8-0623c moved to (I, b) = (25?57, 0?32). 

ses H.E.S.S. measured an extension of a = 0?12 ± 0?02 and 0?05 ± 0?02 of the coincident 

869 TeV source HESS J1837— 069 along the semi-major and semi-minor axis when fitting this 

870 source with an elliptical 2D Gaussian surface brightness profile. This corresponds to a 68% 

871 containment radius of r 68 = 0?18 ± 0?03 and 0?08 ± 0?03 along the semi- major and semi- 

872 minor axis. The size is not significantly different from the LAT 68% containment radius of 

873 r 68 = 0?27±0?07 (less than 2a). Figure [20] shows that the spectrum of 2FGL J1837.3— 0700c 

874 at GeV energies connects to the spectrum of HESS J1837— 069 at TeV energies. 
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888 



HESS J1837— 69 is coincident with the hard and steady X-ray source AX J1838. 0—0655 
( IBamba et al.l 120031 ) . This source was discovered by RXTE to be a pulsar (PSR J1838- 
0655) sufficiently luminous to power the TeV emission and was resolved by Cha ndra to be a 
bright point-like source surrounded by a ~ 2' nebula (IGotthelf fc Halpernll2008l ). The 7-ray 
emission may be powered by this pulsar. The hard spectral index and spatial extension of 
2FGL J1837.3— 0700c disfavor a pulsar origin of the LAT emission and suggest instead that 
the GeV and TeV emission both originate from the pulsa r's wind. There is another X-ray 



point-like source AX J1837.3-0652 near HESS J 18 37-069 (IBamba et a 



resolved into a point-like and diffuse component (IGotthelf fc Halpern 



2003J) that was also 
2008h . Although no 



pulsations have been detected from it, it could also be a pulsar powering some of the 7-ray 
emission. 



10.9. 2FGLJ2021. 5+4026 



887 The source 2FGL J2021. 5+4026 is associated with the 7-Cygni SNR and has been specu- 

888 lated to ori ginate from the interaction of acc elerated particles in the SNR with dense molecu- 

889 lar clouds (|Pollocklll985l ; iGaisser et al.lll998l ). This association was disfa yored when the GeV 

890 emission from this source was detected to be pulsed (PSR J 2021+4026 jAbdo et al.ll2010kl ). 

891 This pulsar was also observed by AGILE (IChen et al.l 120 111 ). 



892 Looking at the same region at energies above 10 GeV, the pulsar is no longer significant 

893 but we instead found in our search an extended source candidate. Figure [25] shows a counts 

894 map of t his source and over lays radio contours of 7-Cygni from the Canadian Galactic Plane 

895 Survey (ITaylor et alj 120031 ) . There is good spatial overlap between the SNR and the GeV 



896 emission. 



897 There is a nearby source 2FGL J2019. 1+4040 that correlates with the radio emission 

898 of 7-Cygni and at GeV energies does not appear to represent a physically distinct source. 

899 We concluded that it was included in the 2FGL catalog to compensate for residuals induced 

900 by not modeling the extension of 7-Cygni and removed it from our model of the sky. With 

901 this model, 2FGL J2021. 5+4026 was found to have an extension a = 0?63 ± 0?05 stat ± 0?04 sys 

902 with TScxt = 128.9. Figure [T71 shows its spectrum. The inferred size of this source at GeV 

903 energies well matches the radi o size of 7-Cygni. M ilagro detected a 4.2a excess at energies 

904 ~ 30 TeV from this location ( Abdo et al. 2009bl lcl). VERITAS also detected an extended 
90s source VER J2019+407 coincident with the SNR above 2 00 GeV and sug gested that the TeV 
906 emission could be a shock-cloud interaction in 7-Cygni ( jWeinsteinl 120091 ) . 
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w 11. Discussion 

9ob Twelve extended sources were included in the 2FGL catalog and two additional extended 

909 sources were studied in dedicated publications. Using 2 years of LAT data and a new analysis 

910 method, we presented the detection of seven additional extended sources. We also reanalyzed 

911 the spatial extents of the twelve extended sources in the 2FGL catalog and the two additional 

912 sources. The 21 extended LAT sources are located primarily along the Galactic plane and 

913 their locations are shown in Figure [26j Most of the LAT-detected extended sources are 

914 expected to be of Galactic origin as the distances of extragalactic sources (with the exception 

915 of the local group Galaxies) are typically too large to be able to resolve them at 7-ray energies. 

916 For the LAT extended sources also seen at TeV energies, Figure [27| shows that there is 

917 a good correlation between the sizes of the sources at GeV and TeV energies. Even so, the 

918 sizes of PWNe are expected to vary across the GeV and TeV ener gy range and the siz e of 



919 HESS J1825— 137 is significantly larger at GeV than TeV energies ( iGrondin et al.ll201ll ). It 

920 is interesting to compare the sizes of other PWN candidates at GeV and TeV energies, but 

921 definitively measuring a difference in size would require a more in-depth analysis of the LAT 

922 data using the same elliptical Gaussian spatial model. 

923 Figure [2S] compares the sizes of the 21 extended LAT sources to the 42 extended H.E.S.S. 



924 sourceso Because of the large field of view and all-sky coverage, the LAT can more easily 

925 measure larger sources. On the other hand, the better angular resolution of air Cherenkov 

926 detectors allows them to measure a population of extended sources below the resolution limit 

927 of the LAT (currently about ~ 0?2). Fermi has a 5 year nominal mission lifetime with a 
92B goal of 10 years of operation. As Figure [TU1 shows, the low background of the LAT at high 

929 energies allows its sensitivity to these smaller sources to improve by a factor greater than the 

930 square root of the relative exposures. With increasing exposure, the LAT will likely begin 

931 to detect and resolve some of these smaller TeV sources. 

932 Figure |2H] compares the spectral indices of LAT detected extended sources and of all 

933 sources in the 2FGL catalog. This, and Tables [3] and HI show that the LAT observes a popula- 

934 tion of hard extended sources at energies above 10 GeV. Figure [201 shows that the spectra of 

935 four of these sources (2FGL J1615.0— 5051, 2FGL J1615. 2-5138, 2FGL J1632.4-4753c, and 

936 2FGL J1837.3— 0700c) at GeV energies connects to the spectra of t heir H.E.S.S. counter parts 

937 at TeV energies. This is also true o f Vela Jr., HESS J1825— 137 fjGrondin et al.ll201ll ). and 



93B RX J1713.7-3946 flAbdo et allEOllt ). It is likely that the GeV and TeV emission from these 



5 The TeV extension of the 42 extended H.E.S.S. sources comes from the H.E.S.S. Source Catalog 



http : / / www . mpi-hd . mpg . de/hf m/HESS/pages/home/sources/. 
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939 sources originates from the same population of high-energy particles. 

940 Many of the TeV-detected extended sources now seen at GeV energies are currently 

941 unidentified and further multiwavelength follow-up observations will be necessary to un- 

942 derstand these particle accelerators. Extending the spectra of these TeV sources towards 

943 lower energies with LAT observations may help to determine the origin and nature of the 

944 high-energy emission. 
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Fig. 1.— Counts maps and TS profiles for the SNR IC 443. (a) TS vs. extension of 
the source, (b) TS ex t for individual energy bands, (c) observed radial profile of counts in 
comparison to the expected profiles for a spatially extended source (solid and colored red 
in the online version) and for a point-like source (dashed and colored blue in the online 
version), (d) smoothed counts map after subtraction of the diffuse emission compared to the 
smoothed LAT PSF (inset). Both were smoothed by a 0?1 2D Gaussian kernel. Plots (a), 
(c), and (d) use only photons with energies between 1 GeV and 100 GeV. Plots (c) and (d) 
include only photo ns which converted in the front part of the tracker and have an improved 
angular resolution (lAtwood et al.ll2009a ). 
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Fig. 2.— A comparison of a 2D Gaussian and uniform disk spatial model of extended 
sources before and after convolving with the PSF for two energy ranges. The solid black line 
is the PSF that would be observed for a power-law source of spectral index 2. The dashed 
line and the dash-dotted lines are the brightness profile of a Gaussian with r68 = 0?5 and the 
convolution of this profile with the LAT PSF respectively (colored red in the online version). 
The dash-dot-dotted and the dot-dotted lines are the brightness profile of a uniform disk 
with r68 = 0?5 and the convolution of this profile with the LAT PSF respectively (colored 
blue in the online version). 
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Table 1. Monte Carlo Spectral Parameters 



Spectral Index 


FluxW 


^l-lOOGeV 


(TS)i_iooGcV ^10- 


lOOGeV 


(TS)lO-lOOGeV 




(ph cm -2 s — 1 ) 










Isotropic Background 


1.5 


3 X IO" 7 
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22233 


18938 


8084 




io- 7 


19079 


5827 


19079 


2258 




3 x icr 8 
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1276 


19303 


541 




IO" 8 


19385 


303 


19381 


112 




3 x 10~ 9 


18694 


62 


12442 


43 
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icr e 


18760 


22101 


18760 


3033 




3 x 10~ 7 
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4913 


18775 


730 




IO" 7 


18804 


1170 


18803 


192 




3 x icr 8 


18836 


224 


15256 


50 




IO" 8 


17060 


50 






2.5 


3 x icr e 


18597 


19036 


18597 


786 




IO" 6 


18609 


4738 


18608 


208 




3 x icr 7 


18613 
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15958 


53 




io- 7 


18658 


203 








3 x IO" 8 


14072 


41 
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icr 5 
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19466 


18354 


215 




3 x 10" 6 


18381 


4205 


15973 


54 




IO" 6 


18449 


966 








3 x icr 7 


18517 
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IO" 7 


13714 


11 






Galactic Diffuse and Isotropic Background' 6 ) 


1.5 


2.3 x icr 8 


90741 


63 






2 


1.2 x 10" 7 


92161 


60 






2.5 


4.5 x 10" 7 


86226 


17 






3 


2.0 x 10~ 6 


94412 


61 







M Integral 100 MeV to 100 GeV flux. 

( b )For the Galactic simulations, the quoted fluxes are the fluxes for sources placed in the Galactic 
center. The actual fluxes are scaled by Equation [12] 



Note. — A list of the spectral models of the simulated point-like sources which were tested for 
extension. For each model, the number of statistically independent simulations and the average value 
of TS is also tabulated. The top rows are the simulations on top of an isotropic background and the 
bottom rows are the simulations on top of the Galactic diffuse and isotropic background. 



-41 - 




Fig. 3. — Cumulative distribution of the TS for the extension test when fitting simulated 
point-like sources in the 1 GeV to 100 GeV energy range. The four plots represent simulated 
sources of different spectral indices and the different lines (colored in the online version) 
represent point-like sources with different 100 MeV to 100 GeV integral fluxes. The dashed 
line (colored red) is the cumulative density function of Equation [TTJ 



-42 - 




Fig. 4.- 



The same plot as Figure [3] but fitting in the 10 GeV to 100 GeV energy range. 
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Fig. 5. — Cumulative distribution of TS cxt for sources simulated on top of the Galactic 
diffuse and isotropic background. 
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Fig. 6. — The distribution of TS values when fitting 985 statistically independent simulations 
of W44. (a) is the distribution of TS values when fitting W44 as a point-like source and (b) is 
the distribution of TS cxt when fitting the source with a uniform disk or a radially-symmetric 
Gaussian spatial model, (c) is the distribution of the change in TS when fitting the source 
with an elliptical disk spatial model compared to fitting it with a radially-symmetric disk 
spatial model and (d) when fitting the source with an elliptical ring spatial model compared 
to an elliptical disk spatial model. 
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Fig. 7. — The distribution of fit parameters for the Monte Carlo simulations of W44. The 
plots show the distribution of best fit (a) flux (b) spectral index and (c) 68% containment 
radius. The dashed vertical lines represent the simulated values of the parameters. 
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Fig. 8. — The detection threshold to resolve an extended source with a uniform disk model 
for a two-year exposure. All sources have an assumed power-law spectrum and the different 
line styles (colors in the electronic version) correspond to different simulated spectral indices. 
The lines with no markers correspond to the detection threshold using photons with energies 
between 100 MeV and 100 GeV, while the lines with star-shaped markers correspond to the 
threshold using photons with energies between 1 GeV and 100 GeV. 
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Fig. 9. — The LAT detection threshold for four spectral indices and three backgrounds (lx, 
10 x, and 100 x the Sreekumar-like isotropic background) for a two-year exposure. The left- 
hand plots are the detection threshold when using photons with energies between 1 GeV 
and 100 GeV and the right-hand plots are the detection threshold when using photons with 
energies between 10 GeV and 100 GeV. The flux is integrated only in the selected energy 
range. Overlaid on this plot are the LAT-detected extended sources placed by the magnitude 
of the nearby Galactic diffuse emission and the energy range they were analyzed with. The 
star-shaped markers (colored red in the electronic version) are sources with a spectral index 
closer to 1.5, the triangular markers (colored blue) an index closer to 2, and the circular 
markers (colored green) an index closer to 2.5. The triangular marker in plot (d) below the 
sensitivity line is MSH 15—52. 



Table 2. Extension Detection Threshold 
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BG 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


E>1 GeV 


1.5 


lx 


148.1 


23.3 


11.3 


8.0 


7.2 


6.9 


6.7 


6.8 


7.1 


7.4 


7.6 


7.9 


8.1 


8.5 


9.2 


9.9 


9.1 


9.2 


9.0 


10.3 




lOx 


148.4 


29.0 


18.7 


15.2 


15.4 


15.0 


16.1 


16.0 


16.8 


17.7 


18.2 


19.3 


20.9 


22.5 


23.8 


24.8 


21.3 


22.8 


23.4 


23.7 




lOOx 


186.8 


55.0 


43.4 


40.7 


41.0 


41.8 


40.9 


40.9 


42.7 


43.6 


38.4 


39.9 


40.6 


38.4 


36.9 


36.3 


37.1 


38.8 


37.2 


37.6 
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lx 


328.4 


43.4 


18.9 


13.4 


11.2 


10.4 


10.2 


10.2 


10.2 


10.4 


10.7 


10.9 


11.2 


11.5 


12.4 


12.6 


13.0 


13.4 


14.0 


14.4 




lOx 


341.0 


55.9 


32.3 


27.6 


26.5 


25.4 


25.6 


25.9 


27.4 


26.8 


27.8 


28.7 


29.8 


30.1 


31.0 


31.5 


31.7 


34.0 


34.3 


35.9 




lOOx 


420.5 


128.3 


90.2 


77.3 


73.3 


70.8 


67.5 


64.3 


64.2 


64.1 


62.8 


63.6 


61.7 


61.9 


58.4 


59.0 


61.4 


63.3 


60.1 


58.1 


2.5 


lx 


627.1 


75.6 


29.8 


19.3 


15.5 


13.5 


12.8 


12.6 


12.5 


12.5 


12.6 


12.9 


12.9 


13.1 


13.5 


13.7 


14.3 


14.8 


15.2 


15.8 




lOx 


638.9 


99.1 


52.1 


39.1 


34.6 


33.0 


32.5 


32.5 


32.8 


33.2 


34.1 


34.3 


34.5 


35.1 


36.6 


36.9 


35.5 


36.0 


36.5 


37.3 




lOOx 


795.0 


262.1 


140.9 


104.3 


90.4 


81.2 


77.2 


75.1 


69.7 


70.9 


66.5 


65.6 


64.9 


64.0 


58.9 


58.1 


60.2 


58.4 


57.5 


55.8 


3 


lx 


841.5 


110.6 


43.2 


25.5 


18.7 


16.1 


14.4 


13.6 


13.3 


13.2 


13.1 


13.1 


13.4 


13.6 


13.5 


13.8 


14.2 


14.4 


14.8 


15.4 




lOx 


921.6 


151.3 


69.1 


47.8 


40.7 


37.1 


35.5 


34.5 


35.1 


35.5 


35.3 


35.3 


35.4 


35.5 


36.8 


37.6 


35.3 


35.4 


36.3 


36.6 




lOOx 


1124.1 


282.9 


181.1 


119.8 


100.7 


91.1 


84.3 


77.9 


73.3 


71.8 


67.6 


66.4 


65.5 


63.9 


59.0 


58.6 


58.8 


57.5 


55.4 


54.4 



E>10 GeV 



lx 


44.6 


8.0 


1.3 


3.2 


2.7 


2.6 


2.5 


2.5 


2.1 


2.5 


2.5 


2.6 


2.7 


2.8 


2.9 


2.9 


3.1 


3.2 


3.3 


3.4 


lOx 


45.2 


9.2 


5.8 


5.0 


4.9 


4.9 


5.0 


5.2 


5.3 


5.7 


5.9 


6.3 


6.6 


6.5 


6.8 


7.6 


7.8 


8.2 


8.5 


8.7 


lOOx 


47.3 


13.4 


11.6 


10.6 


10.8 


10.8 


12.0 


12.7 


13.2 


13.7 


15.3 


16.1 


17.2 


18.2 


18.9 


19.5 


20.4 


21.0 


21.7 


22.9 


lx 


49.7 


8.4 


4.4 


3.3 


2.8 


2.6 


2.6 


2.6 


2.6 


2.6 


2.7 


2.7 


2.8 


2.9 


3.0 


3.2 


3.2 


3.4 


3.5 


3.5 


lOx 


48.6 


9.5 


6.0 


5.2 


5.0 


5.2 


5.2 


5.3 


5.4 


5.8 


6.4 


6.6 


7.0 


7.1 


7.5 


8.0 


8.3 


8.6 


9.0 


9.2 


lOOx 


51.8 


14.7 


11.8 


11.5 


11.5 


11.9 


13.2 


14.0 


14.3 


15.3 


16.2 


16.9 


18.4 


19.2 


19.8 


21.0 


22.0 


22.8 


23.2 


24.3 


lx 


53.1 


9.1 


4.5 


3.3 


2.8 


2.7 


2.6 


2.5 


2.5 


2.6 


2.7 


2.7 


2.8 


2.8 


2.9 


3.1 


3.2 


3.3 


3.5 


3.6 


lOx 


53.7 


10.5 


6.3 


5.4 


5.1 


5.1 


5.3 


5.4 


5.7 


6.0 


6.3 


6.6 


6.8 


6.9 


7.5 


8.1 


8.3 


8.6 


8.9 


9.2 


lOOx 


57.0 


15.6 


12.7 


11.9 


11.8 


12.2 


13.1 


14.3 


14.6 


15.2 


16.3 


17.0 


18.8 


19.2 


19.9 


21.0 


21.9 


22.3 


23.3 


23.7 


lx 


55.5 


9.4 


4.8 


3.4 


2.9 


2.7 


2.6 


2.5 


2.5 


2.5 


2.6 


2.7 


2.7 


2.8 


2.9 


3.0 


3.1 


3.2 


3.4 


3.4 


lOx 


56.0 


10.5 


6.2 


5.3 


5.1 


5.1 


5.1 


5.3 


5.5 


5.7 


5.9 


6.4 


6.4 


6.6 


7.0 


7.8 


8.0 


8.3 


8.6 


8.9 


lOOx 


60.3 


16.2 


12.7 


11.7 


11.8 


12.2 


12.6 


13.8 


14.2 


14.6 


15.8 


16.5 


17.6 


18.5 


19.4 


19.8 


20.7 


21.0 


21.8 


22.5 



Note. — The detection threshold to resolve spatially extended sources with a uniform disk spatial model for a two-year exposure The threshold is calculated for 
sources of varying energy ranges, spectral indices, and background levels. The sensitivity was calculated against a Sreekumar-like isotropic background and the second 
column is the factor that the simulated background was scaled by. The remaining columns are varying sizes of the source. The table quotes integral fluxes in the 
analyzed energy range (1 GeV to 100 GeV or 10 GeV to 100 GeV) in units of 10~ 10 phcm - 2 s~ 1 . 
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Fig. 10. — The projected detection threshold of the LAT to extension after 10 years for a 
power-law source of spectral index 2 against 10 times the isotropic background in the energy 
range from 1 GeV to 100 GeV (solid line colored red in the electronic version) and 10 GeV 
to 100 GeV (dashed line colored blue). The shaded gray regions represent the detection 
threshold assuming the sensitivity improves from 2 to 10 years by the square root of the 
exposure (top edge) and linearly with exposure (bottom edge). The lower plot shows the 
factor increase in sensitivity. For small extended sources, the detection threshold of the 
LAT to the extension of a source will improve by a factor larger than the square root of the 
exposure. 
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Fig. 11. — (a) and (b) are the distribution of TS ext and of TS 2p t s when fitting simulated 
spatially extended sources of varying sizes as both an extended source and as two point-like 
sources, (c) and (d) are the distribution of TS ex t — TS2 P t s for the same simulated sources, 
(a) and (c) represent sources fit in the 1 GeV to 100 GeV energy range and (b) and (d) 
represent sources fit in the 10 GeV to 100 GeV energy range. In (c) and (d), the plus-shaped 
markers (colored red in the electronic version) are fits where TS ext > 16. 
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Fig. 12. — The distribution of TS ext — TS 2p t s when fitting two simulated point-like sources 
of varying separations as both an extended source and as two point-like sources, (a), and 
(b) represent simulations of two point-like sources with the same spectral index and (c) and 
(d) represent simulations of two point-like sources with different spectral indices, (a) and (c) 
fit the simulated sources in the 1 GeV to 100 GeV energy range and (b) and (d) fit in the 
10 GeV to 100 GeV energy range. The plus-shaped markers (colored red in the electronic 
version) are fits where TS cxt > 16. 



-52 - 




Fig. 13. — The cumulative density of TS ex t for the 733 clean AGN in 2LAC that were 
significant above 1 GeV calculated with point like (dashed line colored blue in the electronic 
version) and with gtlike (solid line colored black). AGN are too far and too small to be 
resolved by the LAT. Therefore, the cumulative density of TS cxt is expected to follow a x?/2 
distribution (Equation [11} the dash-dotted line colored red). 



Table 3. Analysis of the twelve extended sources included in the 2FGL catalog 



Name 


GLON 


GLAT 


a 


TS 


TS cxt 


Pos Err 


F1ik("» 


Index 




(cleg.) 


(deg.) 


(deg.) 






(deg.) 






E>1 GeV 


SMC 


302.59 


-44.42 


1.32 ±0.15 ±0.31 


95.0 


52.9 


0.14 


2.7 ±0.3 


2.48 ±0.19 


LMC 


279.26 


-32.31 


1.37 ± 0.04 ±0.11 


1127.9 


909.9 


0.04 


13.6 ±0.6 


2.43 ±0.06 


IC 443 


189.05 


3.04 


0.35 ±0.01 ±0.04 


10692.9 


554.4 


0.01 


62.4 ± 1.1 


2.22 ±0.02 


Vela X 


263.34 


-3.11 


0.88 












Centaurus A 


309.52 


19.42 


- 10 












W28 


6.50 


-0.27 


0.42 ±0.02 ±0.05 


1330.8 


163.8 


0.01 


56.5 ± 1.8 


2.60 ±0.03 


W30 


8.61 


-0.20 


0.34 ±0.02 ±0.02 


464.8 


76.0 


0.02 


29.1 ±1.5 


2.56 ±0.05 


W44 


34.69 


-0.39 


0.35 ±0.02 ±0.02 


1917.0 


224.8 


0.01 


71.2 ±0.5 


2.66 ±0.00 


W51C 


49.12 


-0.45 


0.27 ±0.02 ±0.04 


1823.4 


118.9 


0.01 


37.2 ± 1.3 


2.34 ±0.03 


Cygnus Loop 


74.21 


-8.48 


1.71 ±0.05 ±0.06 


357.9 


246.0 


0.06 


11.4 ±0.7 


2.50 ±0.10 


E>10 GeV 


MSH15-52( fe ) 


320.39 


-1.22 


0.21 ±0.04 ±0.04 


76.3 


6.6 


0.03 


0.6 ±0.1 


2.20 ±0.22 


HESSJ1825-137W 


17.56 


-0.47 


0.65 ±0.04 ±0.02 


59.7 


33.8 


0.05 


1.6 ±0.2 


1.63 ±0.22 



( a ) Integral Flux in units of 10 9 phcm 2 s 1 and integrated in the fit energy range (either 1 GeV to 100 GeV or 
10 GeV to 100 GeV). 



( b >The discrepan c y in t he best fit spectra of MSH 15-52 and HESS J1825-137 compared to lAbdo et all (|2010al ) 
and lGrondin et al.l (|201ll ) is due to fitting over a different energy range. 



Note. — All sources were fit using a radially-symmetric uniform disk spatial model. GLON and GLAT are 
Galactic longitude and latitude of the best fit extended source respectively. The first error on a is statistical and 
the second is systematic (see Section [5]). The errors on the integral fluxes and the spectral indices are statistical 



only. Pos Err is the error on the position of the source. Vela X and the Centaurus A Lobes were not fit in our 
analysis but are included for completeness. 
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Fig. 14. — A TS map generated for the region around the SNR IC 443 using photons with 
energies between 1 GeV and 100 GeV. (a) TS map after subtracting IC 443 modeled as a 
point-like source, (b) same as (a), but IC 443 modeled as an extended source. The cross 
represents the best fit position of IC 443. 
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Fig. 15. — A diffuse-emission-subtracted 1 GeV to 100 GeV counts map of the region around 
2FGL J1856.2+0450c smoothed by a 0?1 2D Gaussian kernel. The plus-shaped marker and 
circle (colored red in the online version) represent the center and size of the source fit with a 
radially-symmetric uniform disk spatial model. The black crosses represent the positions of 
other 2FGL sources. The extension is statistically significant, but the extension encompasses 
many 2FGL sources and the emission does not look to be uniform. Although the fit is 
statistically significant, it likely corresponds to residual features of inaccurately modeled 
diffuse emission picked up by the fit. 



Table 4. Extension fit for the nine additional extended sources 



Name 




GLON 


GLAT 


a 


TS 


TS C xt 


Pos Err 


FluxM 


Index 


Counterpart 






(deg.) 


(deg.) 


(deg.) 






(deg.) 










E>1 GeV 


2FGL J0823.0- 


-4246 


260.32 


-3.28 


0.37 ±0.03 ±0.02 


322.2 


48.0 


0.02 


8.4 ±0.6 


2.21 ±0.09 


Puppis A 




2FGL J1627.0- 


-2425c 


353.07 


16.80 


0.42 ±0.05 ±0.16 


139.9 


32.4 


0.04 


6.3 ±0.6 


2.50 ±0.14 


Ophiuchus 




E>10 GeV 


2FGL J0851.7- 


-4635 


266.31 


-1.43 


1.15 ±0.08 ±0.02 


116.6 


86.8 


0.07 


1.3 ±0.2 


1.74 ±0.21 


Vela Jr. 




2FGL J1615.0- 


-5051 


332.37 


-0.13 


0.32 ±0.04 ±0.01 


50.4 


16.7 


0.04 


1.0 ±0.2 


2.19 ±0.28 


HESS J1616- 


-508 


2FGL J1615.2- 


-5138 


331.66 


-0.66 


0.42 ±0.04 ±0.02 


76.1 


46.5 


0.04 


1.1 ±0.2 


1.79 ±0.26 


HESS J1614- 


-518 


2FGL J1632.4- 


-4753c 


336.52 


0.12 


0.35 ±0.04 ±0.02 


64.4 


26.9 


0.04 


1.4 ±0.2 


2.66 ± 0.30 


HESS J1632- 


-478 


2FGL J1712.4- 


-3941^ 


347.26 


-0.53 


0.56 ±0.04 ±0.02 


59.4 


38.5 


0.05 


1.2 ±0.2 


1.87 ±0.22 


RXJ1713.7- 


3946 


2FGLJ1837.3- 


-0700c 


25.08 


0.13 


0.33 ±0.07 ±0.05 


47.0 


18.5 


0.07 


1.0 ±0.2 


1.65 ±0.29 


HESS J1837- 


-069 


2FGLJ2021. 5+4026 


78.24 


2.20 


0.63 ±0.05 ±0.04 


237.2 


128.9 


0.05 


2.0 ±0.2 


2.42 ±0.19 


7-Cygni 





^Integral Flux in units of 10 9 phcm 2 s 1 and integrated in the fit energy range (either 1 GeV to 100 GeV or 10 GeV to 100 
GeV). 



(^The discrepancy in the best fit spectra of 2FGL J1712.4— 3941 compared to lAbdo et al.l ( 201ll ) is due to fitting over a different 
energy range. 



Note. — The columns in this t able have the same meaning as those in Table RX J1713.7— 3946 and Vela Jr. were previously 



studied in dedicated publications (jAbdo et al. 



2011 



Tanaka et al. 



2011 



Table 5. Dual localization, alternative PSF, and alternative approach to modeling the diffuse emission 



Name 


TSp i n tlike 


TSgtiike 


TS a i t) diff 


TS cx t point like 


TSextgtlike 


TS ex talt,diff 


a 


falt.diff 


falt,psf 


TS2pts 
















(deg.) 


(deg.) 


(deg.) 




E>1 GeV 


2FGL J0823.0-4246 


331.9 


322.2 


356.0 


60.0 


48.0 


56.0 


0.37 


0.39 


0.39 


23.0 


2FGL J1627.0— 2425c 


154.8 


139.9 


105.7 


39.4 


32.4 


24.8 


0.42 


0.40 


0.58 


24.5 


E>10 GeV 


2FGL J0851. 7-4635 


115.2 


116.6 


123.1 


83.9 


86.8 


89.8 


1.15 


1.16 


1.17 


15.5 


2FGL J1615.0— 505l( a ) 


48.2 


50.4 


56.6 


15.2 


16.7 


17.8 


0.32 


0.33 


0.32 


13.1 


2FGLJ1615. 2-5138 


75.0 


76.1 


83.8 


42.9 


46.5 


54.1 


0.42 


0.43 


0.43 


35.1 


2FGL J1632.4-4753c 


64.5 


64.4 


66.8 


23.0 


26.9 


25.5 


0.35 


0.36 


0.37 


10.9 


2FGL J1712. 4-3941 


59.8 


59.4 


39.9 


38.4 


38.5 


30.7 


0.56 


0.55 


0.53 


2.7 


2FGL J1837.3-0700c 


44.5 


47.0 


39.2 


17.6 


18.5 


16.1 


0.33 


0.32 


0.38 


10.8 


2FGLJ2021. 5+4026 


239.1 


237.2 


255.8 


139.1 


128.9 


138.0 


0.63 


0.65 


0.59 


37.3 



W Using pointlike, TS ex t for 2FGL J1615. 0—5051 was sligthly below 16 when the source was fit in the 10 GeV to 100 GeV energy 
range. To confirm the extension measure, the extension was refit in pointlike using a slightly lower energy. In the 5.6 GeV to 100 
GeV energy range, we obtained a consistent extension and TS cxt =28.0. In the rest of this paper, we quote the E > lOGeV results for 
consistency with the other sources. 

Note. — TSp i n tiikei TSgtiike > and TS a it,difT are the test statistic values from pointlike, gtlike, and gtlike with the alternative 
approach to modeling the diffuse emission respectively. TS ex tpointiike > TS cx t g tiike> an d TS ex t a it,diff are the TS values from pointlike, 



gtlike, and gtlike with the alternative approach to modeling the diffuse emission respectively, a, c" a it,difF? an d & a \ t ps f are the fit sizes 
assuming a radially-symmetric uniform disk model with the standard analysis, the alternative approach to modeling the diffuse emission, 
and the alternative PSF respectively. 
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Table 6. Nearby Residual-induced Sources 



Extended Source 



Residual-induced Sources 



2FGL J0823. 0—4246 
2FGL J1627.0— 2425c 
2FGLJ0851. 7-4635 
2FGL J1615. 0—5051 
2FGL J1615.2— 5138 
2FGL J1632.4— 4753c 
2FGL J1712.4— 3941 
2FGL J1837.3— 0700c 
2FGL J2021.5+4026 



2FGL J0821. 0-4254, 2FGL J0823.4-4305 

2FGL J0848.5— 4535, 2FGL J0853.5— 4711, 2FGL J0855.4— 4625 

2FGL J1614.9— 5212 
2FGL J1634.4-4743c 

2FGL J1835.5— 0649 
2FGL J2019. 1+4040 



Note. - For each new extended source, we list nearby 2FGL soruces that we have 
concluded here correspond to residuals induced by not modeling the extensions of nearby 
extended sources. 



-61 - 




Fig. 16. — A diffuse-emission-subtracted 1 GeV to 100 GeV counts map of 
2FGL J0823. 0—4246 smoothed by a 0?1 2D Gaussian kernel. The triangular marker (col- 
ored red in the online version) represents the 2FGL position of this source. The plus-shaped 
marker and the circle (colored red) represent the best fit position and extension of this source 
assuming a radially-symmetric uniform disk model. The two star-shaped markers (colored 
green) represent 2FGL sources that were removed from the background model. From left 
to right, these sources are 2FGL J0823.4— 4305 and 2FGL J0821. 0—4254. The lower right 
inset is the model predicted emission from a point-like source with the same spectrum as 
2FGL J0823.4— 4305 smoothed by the same kernel. This source is spatially coincident with 
the Puppis A SNR. The light blue co ntours correspond to the X-ray image of Puppis A 
observed by ROSAT jPetre et al.lll99fih . 
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Fig. 17.— The spectral energy distribution of the extended sources Puppis A 
(2FGL J0823. 0—4246) and 7-Cygni (2FGL J2021. 5+4026). The lines (colored red in the 
online version) are the best fit power-law spectral models of these sources. Puppis A has a 
spectral index of 2.21 ± 0.09 and 7-Cygni has an index of 2.42 ± 0.19. The spectral errors 
are statistical only. The upper limit is at the 95% confidence level. 
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Fig. 18. — A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of 
2FGL J0851.7— 4635 smoothed by a 0?25 2D Gaussian kernel. The triangular marker (col- 
ored red in the electronic version) represents the 2FGL position of this source. The plus- 
shaped marker and the circle (colored red) are the best fit position and extension of this 
source assuming a radially-symmetric uniform disk model. The three black crosses repre- 
sent background 2FGL sources. The three star-shaped markers (colored green) represent 
other 2FGL sources that were removed from the background model. They are (from left 
to right) 2FGL J0853.5— 4711, 2FGL J0855.4— 4625, and 2FGL J0848.5— 4535. The circular 
and square-shaped marker (colored blue) represents the 2FGL and relocalized position of 
another 2FGL source. This extended source is spatially coincident with t he Vela Jr. SNR. 



The co ntours (colored light blue) correspond to the TeV image of Vela Jr. ( lAharonian et al. 
2007ah . 
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Fig. 19.— A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of 
2FGL J1615. 0—5051 (upper left) and 2FGL J1615.2— 5138 (lower right) smoothed by a 0?1 
2D Gaussian kernel. The triangular markers (colored red in the electronic version) represent 
the 2FGL positions of these sources. The cross-shaped markers and the circles (colored red) 
represent the best fit positions and extensions of these sources assuming a radially sym- 
metric uniform disk model. The two black crosses represent background 2FGL sources and 
the star-shaped marker (colored green) represents 2FGL J1614. 9-5212, another 2FGL source 
that was removed from the background model. The contours (colored light blue) correspond 



to th e TeV image of HESS J1616-508 (left) and HESS J1614-518 (right) flAharonian et al. 
2006h . 
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Fig. 20.— The spectral energy distribution of four extended sources associated with 
unidentified extended TeV sources. The black points with circular markers are obtained 
by the LAT. The points with plus-shaped markers (colored red in the electronic ver- 
sion) are for the associated H.E.S.S sources, (a) the LAT SED of 2FGL J1615. 0—5051 
together with the H.E.S.S. SED of HESS J1616-508. (b) 2FGL J1615.2— 5138 and 
HESS J1614-518. (c) 2FGL J1632.4-4753c and HESS J1 632-478. fd) 2FGL J18 37.3-0700c 
and HESS J1837-069. The H.E.S.S. data points are from JAharonian et alihoOfil ). Both LAT 
and H.E.S.S. spectral errors are statistical only. 
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Fig. 21. — A diffuse-emission-subtracted 1 GeV to 100 GeV counts map of (a) the region 
around 2FGL J1627. 0—2425 smoothed by a 0?1 2D Gaussian kernel and (b) with the emission 
from 2FGL J1625.7— 2526 subtracted. The triangular marker (colored red in the online 
version) represents the 2FGL position of this source. The plus-shaped marker and the circle 
(colored red) represent the best fit position and extension of this source assuming a radially- 
symmetric uniform disk model and the black cross represents a backgr ound 2FGL source . 
The contours in (a) correspond to the 100 jum image observed by IRAS ( lYoung et al.lll986l ). 
The contours in (b) corres pond to CO {J = 1 — > 0) emission integrated from —8 kms" 1 to 
20 kms" 1 . They are f rom Ide Geus et al.l ( 119901 ) . were cleaned using the moment-masking 
technique (jDamell201ll ). and have been smoothed by a 0?25 2D Gaussian kernel. 
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Fig. 22.— A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of 
2FGL J1632.4— 4753c (a) smoothed by a 0?1 2D Gaussian kernel and (b) with the emission 
from the background sources subtracted. The triangular marker (colored red in the elec- 
tronic version) represents the 2FGL position of this source. The plus-shaped marker and the 
circle (colored red) are the best fit position and extension of 2FGL J1632.4— 4753c assuming 
a radially-symmetric uniform disk model. The four black crosses represent background 2FGL 
sources subtracted in (b). The circular and square-shaped markers (colored blue) represent 
the 2FGL and relocalized positions respectively of two additional background 2FGL sources 
subtracted in (b). The star-shaped marker (colored green) represents 2FGL J1634.4— 4743c, 
another 2FGL source that was removed from the background model. The contours ( colored 
light blue) correspond to the TeV image of HESS J1632— 478 ( lAharonian et al.l 120061 ) . 
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Fig. 23.— A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of 
2FGL J1712.4— 3941 (a) smoothed by a 0?15 2D Gaussian kernel and (b) with the emis- 
sion from the background sources subtracted. This source is spatially coincident with 
RX J1713.7— 3946 and was recently studied in lAbdo et al.l (120 111 ). The triangular marker 
(colored red in the online version) represents the 2FGL position of this source. The plus- 
shaped marker and the circle (colored red) are the best fit position and extension of this 
source assuming a radially symmetric uniform disk model. The two black crosses represent 
background 2FGL sources subtracted in (b ). The contours (colored light blue) correspond 
to the TeV image (lAharonian et al.ll2007bl ). 
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Fig. 24. — A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of the region 
around 2FGL J1837.3— 0700c (a) smoothed by a 0?15 2D Gaussian kernel and (b) with the 
emission from the background sources subtracted. The triangular marker (colored red in the 
online version) represents the 2FGL position of this source. The plus-shaped marker and 
the circle (colored red) represent the best fit position and extension of 2FGL J1837.3— 0700c 
assuming a radially-symmetric uniform disk model. The circular and square-shaped markers 
(colored blue) represent the 2FGL and the relocalized positions respectively of two back- 
ground 2FGL sources subtracted in (b). The star-shaped marker (colored green) repre- 
sents 2FGL J1835.5— 0649, another 2FGL source that was removed from the background 
model. The contours (c olored light blue) correspond to the TeV image of HESS J1837— 069 
(lAharonian et al.ll2006l ). The diamond-shaped marker (colored orange) represents the posi- 
tion of PSR J1838— 0655 an d the hexagonal-shaped m arker (colored purple) represents the 
position AX J1837.3— 0652 jGotthelf fe Halpernlhoosl ). 
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Fig. 25. — A diffuse-emission-subtracted 10 GeV to 100 GeV counts map of the region 
around 2FGL J2021. 5+4026 smoothed by a 0?1 2D Gaussian kernel. The triangular marker 
(colored red in the online version) represents the 2FGL position of this source. The plus- 
shaped marker and the circle (colored red) represent the best fit position and extension of 
2FGL J2021. 5+4026 assuming a radially-symmetric uniform disk model. The star-shaped 
marker (colored green) represents 2FGL J2019. 1+4040, a 2FGL source that was removed 
from the background model. 2FGL J2021. 5+4026 is spatially coincident with the 7-Cygni 
SNR. The contours (colored light blue) co rrespond to the 408 MHz image of 7-Cygni observed 
by the Canadian Galactic Plane Survey (ITaylor et al.ll2003l ). 
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Fig. 26. — The 21 spatially extended sources detected by the LAT at GeV energies with 2 
years of data. The twelve extended sources included in 2FGL are represented by the circular 
markers (colored red in the online version). The nine new extended sources are represented 
by the triangular markers (colored orange). The source positions are overlaid on a 100 MeV 
to 100 GeV Aitoff projection sky map of the LAT data in Galactic coordinates. 
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Fig. 27.— A comparison of the sizes of extended sources detected at both GeV 
and TeV energies. The TeV sizes of W30, 2FGL J18 37.3-0700c, 2FGL J163 2.4-4753c, 
2FGL J1615. 0—5051, and 2FGL J1615.2— 5138 are from kharonian et all J2006h . The TeV 
sizes of MSH 15-52, HESS J 1 825-1 37 Vela X Vela Jr., RX J1713.7-3946 and W28 are 
from lAharonian et all (bo05l l2006bl lc3. |2007allbL bood). T he TeV size of IC 443 is from 



Acciari et all (120091 ) and W51C is from 



Krause et all (120 111 ). The TeV sizes of MSH 15-52, 



HESS J1614-518, HESS J1632-478, and HESS J1837-069 have only been reported with an 
elliptical 2D Gaussian fit and so the plotted sizes are the geo metric mean of the semi-major 
and semi- minor axis. The LAT extension of Vela X is from lAbdo et al.l (|2010fl ). The TeV 
sources were fit assuming a 2D Gaussian surface brightness profile so the plotted GeV and 
TeV extensions were first converted to r68 (see Section |2~4"1) . Because of their large sizes, the 
shape of RX J1713.7— 3946 and Vela Jr. were not directly fit at TeV energies and so are 
not included in this comparison. On the other hand, dedicated publications by the LAT 



collaboration on these s ources showed that their morphologies are consistent f lAbdo et al. 



20 111 ; iTanaka et al.ll201ll ). The LAT extension errors are the statistical and systematic errors 



added in quadrature. 
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Fig. 28.— The distributions of the sizes of 18 extended LAT sources at GeV energies 
(colored blue in the electronic version) and the sizes of the 40 extended H.E.S.S. sources 
at TeV energies (colored red). The H.E.S.S. sources were fit with a 2D Gaussian surface 
brightness profile so the LAT and H.E.S.S. s izes were first converted to r 68 . The GeV size 
of Vela X is taken from lAbdo et al.l (l2010il ). Because of their large sizes, the shape of 
RX J1713.7— 3946 and Vela Jr. were not directly fit at TeV energies and are not included in 
this comparison. Centaurus A is not included because of its large size. 
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Fig. 29. — The distribution of spectral indices of the 1873 2FGL sources (colored red in 
the electronic version) and th e 21 spatiall y exte nded sources (colored blue). The index 
of Centaurus A is taken from iNolan et al.l (120121 ) and the index of Vela X is taken from 



Abdo et al.1 fl2010fl ). 



